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ABSTRACT 


This analytical study focuses on the response of unsymmetrically laminated cylinders 
subjected to thermally-induced preloading effects and compressive axial load. Attention 
is focused on the displacement response and three-dimensional stress state of cylinders 
having {+45I-45I0 Jk , [+45/-45/0J fr , and [Of -45 1 + 45] „ stacking sequences with 
clamped end conditions. The methods used in the analyses involve derivation of the plane 
stress and three-dimensional equilibrium equations and boundary conditions using the 
method of minimum total potential energy with nonlinear strain-displacement relations. 
The plane stress equations and boundary conditions are solved in closed-form for the 
displacements and intralaminar stresses. The three-dimensional equilibrium equations 
are then solved for the interlaminar shear stress x xr using the results of the plane stress 
problem. For the three cylinders analyzed, the radial deformations are observed to be 
larger for the unsymmetrically laminated cylinders, particularly in the boundary layer 
near the ends of the cylinders. With the nonlinear effects included, the boundary layer 
length increases with increasingly compressive axial load. If the thermally-induced 
preloading effects are not included, the deformations and intralaminar stresses are under- 
predicted. Also, it is observed that the boundary conditions for the axial load must 
include the thermally deformed shape of the cylinder. At low axial load levels, it was 
seen that both the fiber-direction intralaminar stress and the interlaminar stress x xr are 
dominated by the thermally-induced preloading effects. However , the intralaminar stress 
perpendicular to the fiber direction and the intralaminar shear stress are largely 
unaffected by the thermally-induced preloading effects. 
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I. INTRODUCTION 


Cylinders made of composite materials are structurally efficient and well suited to 
manufacture by automated fiber-placement techniques such as filament winding. 
Composite cylinders are also more tolerant of unsymmetric stacking sequences than are 
flat plates. Unsymmetrically laminated plates warp when cooled from the consolidation 
temperature during manufacture. This effect becomes problematic for the designer, since 
the degree of warp must be anticipated and accounted for before manufacture. Because 
of the difficulties brought about by the manufacture of unsymmetrically laminated plates, 
there is a tendency to view the idea of unsymmetrically laminated cylinders with the 
same skepticism. However, due to the axisymmetric geometry inherent in the cylindrical 
form, cylinders resist much of the warping observed in the manufacture of plates. In 
practice, cylinders which are manufactured with the filament winding process are 
sometimes purposely made with unsymmetric stacking sequences because it is easier and 
less costly to complete the winding process without changing the fiber orientation on the 
outer layers of the cylinder. 

One of the benefits of designing a structure with composite materials is the ability to 
tailor the structural properties of a member not only through changing the material type 
or thickness, but also by taking advantage of the couplings between bending, stretching, 
and shearing possible through changes in the stacking sequence of the wall of the 
structure. These couplings may impart a structure with behaviors which are beneficial 
under certain types of loadings, for example, increased buckling resistance to 
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compressive axial loads. For tubular members such as cylinders, these couplings will 
also impart to the structure a nonuniform shape as a result of cooling from consolidation 
temperature to ambient temperature. Although considerable attention to the response of 
composite cylinders has been given by others in the past, the combined effects of 
thermally-induced preloading effects and a compressive axial load on unsymmetrically 
laminated cylinders have previously received little or no attention. In fact, buckling 
analyses of unsymmetric cylinders are available in the literature which completely and 
erroneously ignore thermally-induced deformation (ref. 1). 

The purpose of this work is to investigate the effect of applying a compressive axial 
load to an unsymmetrically laminated cylinder which has already deformed in a 
nonuniform manner due to a temperature change from the consolidation temperature to 
ambient temperature. The work is important because cylinders subjected to a 
compressive axial load are prone to buckle. Analysis of the buckling phenomenon is a 
difficult problem. However, a key component of the buckling analysis is the prebuckling 
analysis. The predicted prebuckling state for an unsymmetrically laminated cylinder 
deformed by thermal effects before loading could be significantly different than the 
prebuckling state for that cylinder predicted by an analysis which ignores thermal effects. 
Though interest here will not extend to a buckling analysis, as it is considerably beyond 
the scope of the present effort, the formulation addressed here must be considered if a 
buckling analysis is to be developed. The analysis developed here is indeed a 
prebuckling analysis. Implicitly, it will be possible to determine if ignoring thermally- 
induced deformations may have an impact on predicting buckling loads. Towards this 
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end, the equilibrium equations and consistent boundary conditions which govern the 
response of thin cylindrical panels under general loadings will be derived in the second 
chapter. In the third chapter, these equations will be simplified due to conditions of 
axisymmetric geometry and response for the case of compressive axial end loading. The 
equations will be solved, to include the thermally-induced preloading effects, for the 
axial, tangential, and radial displacements as a function of the cylinder’s length 
coordinate. These displacements will be calculated for three cylinders with different 
stacking sequences and are presented in graphical form. 

In the fourth chapter, the solutions derived previously for the displacements will be 
used to obtain relations for the intralaminar stress components within each layer of the 
cylinder. Relations between the stress components and the radial coordinate for clamped 
boundary conditions and several axial load levels will be graphically presented for the 
three cylinders. In the fifth chapter, the equilibrium equations and boundary conditions 
in cylindrical coordinates for a three-dimensional stress state will be presented and 
simplified based on an investigation of the magnitudes of the coupling terms in the 
equations. The sixth chapter will present a solution of the simplified three-dimensional 
equilibrium equations and boundary conditions by making use of the stress-strain 
relations of the fourth chapter, which are of closed form. The significant interlaminar 
stress component will be calculated for each of three cylinders and the results will be 
graphically presented. 
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A. Cylinder Nomenclature and Geometry 


The cylinder is referenced to a rectangular coordinate system with the axis of the 
cylinder centerline coincident with the X axis, as shown in Fig. 1. 



For convenience, the origin of the global coordinate system is chosen to be at the 
midspan of the cylinder. Naturally, a cylindrical coordinate system is used for the 
analysis. The cylindrical coordinate system consists of an x axis which is coincident with 
the X axis, 0 , which is measured positive from the +Z axis toward the + Y axis, and r, 
which is measured outward from the X axis. The cylinder has a length L, a mean radius 
R, and wall thickness H. Within the cylinder wall, a 2 coordinate is defined as being 
positive outward, measured from the mean radius R. The axial, tangential, and radial 
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displacement components of a point in the cylinder wall are denoted by u(x,0,r), v(x,0,r), 
and w(x,0,r), respectively. The axial displacement u(x,0,r) is measured positive along 
the +X axis, the tangential displacement v(x,0,r) is measured positive in the +0 
direction, and the radial displacement w(x,0,r) is measured positive outward in the r 
coordinate direction. 

The derivation of the governing equations follows in the next chapter. 
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H. DERIVATION OF THE EQUILIBRIUM EQUATIONS 
FOR THIN CYLINDRICAL PANELS 


A. The Method of Minimum Total Potential Energy 
The total potential energy of a body is given by 

n = u + n^, 


(i) 


where U is the strain energy of the body given by 

u ■ \ 

x.j and being the stress and strain tensors. The potential energy 

“ n w - ///r'l “.*-///.«. *• 


(2) 

n ^ is defined 
( 3 ) 


where B t is the body force distribution acting on the volume V of the body, T t are the 
surface tractions acting over the boundary S of the body, and u t is the displacement 
field. 

For thin cylindrical shells, the strain energy expression of eq. (1) can be simplified 
by assuming a state of plane stress, i.e., by assuming that o rr , t 6r , and x xr are zero. 
As a further simplification, the body force distribution will also be assumed to be zero. 
Considering preloading deformations and including these two restrictions, eq. (1) can be 
expanded to yield 
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n = \ fffy[( 0 x~ a *) e x + (°e”°e) e e + ( T *e _T Jrt)Yxe]^ 
-ff s (T x u + T 0 v + T z w)dS. 


( 4 ) 


For a cylindrical panel with mean radius R, length L, opening angle p , thickness H, 

bounded by the six surfaces x = x=+— , 0 = , 0 = + -§. r = R-— t and 

2 2 2 2 2 
JJ 

r = /?+— , eq. (4) becomes 
2 


X m +— 

2 2 2 


n= | Iff [(°x - °f) c x + (°e - °0) c e + ( T ie -’'*')y *\ rdQdxdr + n w , (5) 


r-R-— x—- 0— i 
2 2 2 


where 


E loaJ = ~ ff s ( T x U + T 9 V + (6) 

is the potential energy of the applied loads. The stress components superscripted with 
a "P" denote preloading effects. These components could be due to imperfect cylinder 
geometry, thermally-induced deformations, or any other influence unrelated to the 
loading. 

The functional in eqs. (5) and (6) represent the total potential energy of a cylindrical 
panel under the assumption of plane stress. This functional will be minimized in a later 
section through the methods of variational calculus. 
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B. Assumptions Related to Thin Shells 

The following are Donnell’s assumptions for the kinematics of deformation : 

u(x,Q,r) = u °(x,B) + zp"(*,6) 
v(x,0,r) = v°(x,0) + zpo(x,0) 
w(x,0,r) = w°(x,0) . 


(7) 


( 8 ) 


In the above, the local thickness coordinate z is given by 

z = r - R, 

the superscript zero denotes displacements of the cylindrical panel’s reference surface 
taken at the mean radius R, and the 0’s are the rotations at the reference surface given 
by 




dw° 

etc 

dw° 

Rdd 


(9) 


Note that the displacements u°, v° , and w° and the rotations p® and Pe are relative 
to a perfectly cylindrical panel before any preloading effects. The pertinent strain- 
displacement relations in cylindrical coordinates are 


&L + 

If dw 

etc 

2{dx 

dv 

+ ^ + 

RdB 

r 

i*: + 

+ 

etc 

r 30 


dw dw 
rdd etc 


( 10 ) 


For thin cylindrical panels, the approximation 

r = R 


(ID 
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can be made with sufficient accuracy. Substituting eqs. (7) and (9) into eq. (10) using 
eq. (11), the strains become 


O O 

®x = C * + ZK X 
o o 

= ®e + * K e 

= Y^0 + Zi& , 


, _£P? + ap; 

* ax Rde ' 

The stresses within the cylindrical panel are given by 

°x ~ Qn( c , “ c i) + Qi2( c e " c e) + Qw(Yxe ~ Yx«) 

= Qn^®* ~~ ®x) + Q22(®e ~~ *o) + Q26 ^Y,o - Yjfl) 
T x6 “ Ql6 ^® x ~ ®x> + Q26^®0 ” ®e) + Q^Yjt ~ Yjj) 1 

where e x , e£ and y^e are the strains due to preloading effects. 
Equation (14) can be rewritten as 

°x = Qn®* + Qi2®e + QwYjtf ” °x 

°e = Qi2®, + Q22®e + Q2«Y x e-oJ 

T x« = Q|6 ®x + ^26® 8 + Q«Yx 6 " T xfl » 


where 


O 

®* = 


®e ~ 


O 

** = 


du° 

+ — B " 2 

etc 

T* 

dv° 

w° 

+ — — -4 

Rde 

R 

dv° 

+ &1 

dx 

Rde 

dK 

o 

dx 

; *e = 


!»•* 


i«n® 


ap£ . 

Rd& ’ 


where 


( 12 ) 


(13) 


(14) 


(15) 
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( 16 ) 


°r = (Qii £ r + Qi2 e £ + Qi«Y*) 

°0 = (@12 C x + Q22 C 0 + 

T x6 = (@16 e x + Q26 C fl + O^Yxo) 


If the preloading effects are due to thermally-induced deformations (from a 
temperature change such as from consolidation temperature to room temperature, for 
example), then 


«; = «A T 

c e =a e AT 

V* = a xa* T > 


(17) 


where a x , a e and are the coefficients of thermal expansion of the material in the 
cylindrical coordinate system. If this is the case, then the stress-strain relation of eq. 
(14) can be written as 

°x = Qn( e x“ a x A7 ) + Qi 2 ( e e" a e A7 ) + Qi 6^fxa~ a x» AT > 

o e = Q 12 (e,-a,A 7 ) + Q 22 (e 6 -a e A 7 ) +Q 26 (Y jrf T a xe A7 ) (18) 

= Q 16 (e,-a,A7) +Q 26 (e e -« e A *) +Q 66 (Y, 0 -«xe A7 ) , 


or 


where 


° x = Qn c x + Qi2 c e + ^167x6 ” °x 

°e = Ql2 6 x + ^22 c e + “ °® 

T xe = Qi6 e x + Qz6 e e + QwY* - 

°J ■ (Qn a x + Qi2 a e + Qi« a x6) Ar 
°e = (Ql2 a * + ^22*6 + Q26 a xe) A7, 

= (Ql6 a x + Q26°e + Q66«xe) Ar - 


(19) 


( 20 ) 


Here the superscript "T" denotes that the preloading effects are thermally-induced. 
The stresses oj, oj, tJ, would be the stresses at a point arising from a temperature 
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change if the composite is fully constrained from any deformation. Since, in general, 


each layer in the laminated cylindrical panel has unique values for Q (j , a x , a 0 , and a i0 , 
there are separate equations, eqs. (19) and (20), relating the strains e x , e e , and y x0 , and 
the temperature change, AT, to the stresses in each layer. 


C. Specification of the Potential Energy Due to External Loads 

In the absence of body forces, the potential energy due to applied tractions is given 
by 

,v **>•]&, (21) 

where T x , T e , and T z are the known applied tractions acting on the surfaces of the 
cylindrical panel and the tractions and displacements u and v are functions of x, 6, and 
z, and w is a function of x and 6. Substituting eq. (7) into the above, 


= -//, [t;(j:,0,z)(b °(x,0)+zp^(x,e)) 

♦ T Q (x,bj{ V 0 (x,Q)+zVl(x,e)) (22) 

+ r z (x,0^)w*(x,0)]&. 

As shown in Fig. 2, considering the tractions acting along the x=+L/2 edge of the 
cylindrical panel, the contribution towards n frif is 




r *(*M v H e M4 e I 
r -( + HH e )K • 


( 23 ) 
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Fig. 2. Tractions acting on the x = +L/2 edge. 


Since u° , v°, w° , p", and p£ are not functions of z, the integration with respect to 
z can be distributed, resulting in 





+ 






(24) 


The integrals with respect to z are the resultant forces and moments acting along the 
x=+L/2 edge and are defined by 


3 




These resultants have the dimension of force or moment per unit circumferential length 
and they are illustrated in Fig. 3. Substituting the definitions of eq. (25) into eq. (24), 
the contribution towards ^load due to the tractions acting along the x= +L/2 edge can 
be written as 


In a similar fashion, the contribution to n fr i l f due to tractions acting along the x=-L/2 
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edge can be written as 


nU.- r /.lK (6) “1-|’ 6 ) tM;(e),,: (-T e ) 07, 

where the resultant forces and moments along the x=-L/2 edge are given by 
N X '(Q) = ; M' x { 6) 

; w„(9)=/J^-|,e^jofe (28) 

<?;<«, - />{-£*)* • 

These resultants have the dimension of force or moment per unit circumferential length. 
The tractions and associated resultants at x=-L/2 are shown in Fig. 4 and Fig. 5. 



* 
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Referring to Fig. 6 and considering the tractions acting along the 0=+/3/2 edge of 



Since u° , v°, w°, p®, and p£ are not functions of z, the integration with respect to 
z can be distributed, resulting in 
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The integrals with respect to z are the resultant forces and moments acting along the 
0= +0/2 edge. They are defined as 


B M 

KM - u T {*> + ^y ; 

<(*) = f j T^c,^y. ; Mfr) = f j r^+|^jzdz (31) 

Qe( x ) = /_* ' 

These resultants have the dimension of force or moment per unit length. The orientation 
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of these stress resultants is shown in Fig. 7. 



The contribution towards n w due to the tractions acting along the 6= +(5/2 edge 
can therefore be written as 

n a..i * (32) 

+fV 0 + (x)v°|x,+-|j +Af e + (x)p^x,+-|j +<?« (x)w°|x,+-| jjdfc . 

In a similar fashion, the contribution to ^lood due to tractions acting along the 0=-/S/2 
edge can be written as 

n ~l.~f ■/^K W "H) +JWW H) (33) 

+Af e '(x)v + Af 6 '(x)p^.x,--|j , 

where the resultant forces and moments along the 0=-0/2 edge are given by 
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(34) 


HUM - ; UiM - fjTjf- 

# +!L \ 

AW - ; W.W */.| ■ T* 


Fig. 8 and Fig. 9 show the tractions and associated resultants acting on the 8 0/2 edge. 
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Finally, the tractions acting on the top and bottom surfaces of the cylindrical panel 
shall be considered. In this work, only loadings normal to these surfaces will be 
considered, i.e., T x and T e will be taken as zero. These are illustrated in Fig. 10. 





Fig. 10. Tractions acting on the top and bottom surfaces. 

With this limitation in the loading on the top and bottom surfaces, the contribution to 
^load due to normal tractions acting on the top surface is 



L 

-fl C. 1 4 .a,f)w^e^ 


(35a) 


The contribution towards ^leod due to normal tractions acting on the bottom surface is 



+ I 


L.L /,*.!| 


(35b) 


The definitions 


and 


9 + (x,0) S r^e, + |j 


(36a) 
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are introduced here to give a more familiar meaning to these loadings, i.e., that of a 
distributed normal pressure loading with dimensions of force per unit area. Substituting 
eqs. (36) into eqs. (35), 


L 




q*(xfi)w°(xfi)RdQdx , 


(37a) 


and i j> 

n^l,.-i=/3/e'.-i ? te9,w “ te6) ' W0 ' lr • (37b) 

2 2 2 

For convenience, and for ease of discussion, commonly occurring loadings will be 
partitioned from the general form of the potential energy due to surface tractions. For 
the case of known axial loads applied at the x=-L/2 and x= +L/2 edges of the 
cylindrical panel, 

v w;(8)» ”(-f .sjft® n;( 6)« . (38a) 

For bending moments applied along the x=-L/2 and x=+L/2 edges, 

- /.i . < 38 «» 

For inplane shearing loads applied along the x=-L/2 and x= +L/2 edges, and along the 
0=-/J/2 and 0=+/3/2 edges of the cylindrical panel, 


(38c) 


♦1 ♦! 
n i^=/*\W)v‘(-o 9 )W0 -f'\N^B)v‘(+k6)IW 

2 L 2 Z 
+ i. t L 


For shearing moments applied along the x=-L/2 and x=+L/2 edges, and along the 
d=-i 3/2 and 6= +/J/2 edges, 


♦2 


♦2 


n — = Oi - 0-1 "»PS(*f •»)*» 

* O-i - O-i <wp^. + {)* • 

For out-of-plane shear loadings applied along the x=-L/2 and x=+L/2 edges, 


(38d) 


+2 


+2 


n ^Oi ■ O-i WOO 9 }™ 6 


For circumferential loadings applied along the 0=-/3/2 and 0=+(3/2 edges, 


+ 1 + £ 

n w=/ I .!r K °( x> "f * 


For bending moments applied along the 0=-/3/ 2 and 6—+fi/2 edges, 


♦2 *2 

* Oi - Oi "• wp^ + |)* • 


For out-of-plane shear loadings applied along the $=-012 and 0=+/3/2 edges, 


(38e) 


(380 


(38g) 
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(38h) 


+— 

n iood = f x .\ <?8 W w 




For normal distributed loadings q (x,0) and q*(x,G) applied to the bottom and top 
surfaces of the cylindrical panel, respectively, 


n 


load 


♦I 




= f * f \ (q-(x,Q) -q\xfi))w°{x,e)RdQdx 


(38i) 


Note that, in this theory, the difference between the normal distributed loadings q '(*,0) 
and q*(x,6) has an effect, not the individual distributed loadings themselves. In other 
words, distributed normal loadings of q ”(x,0) = ~q and q *(*,0) = 0 would have the same 
effect as distributed normal loadings of q (x,0) = 0 and q *(x>6) = q • Hence, since the 
theory utilizes only the difference in the loadings, only the difference 

„ _ „♦ (39) 

shall be used in the following. The sense of this net distributed load q is outward 
positive and it acts normal to the surfaces of the cylindrical panel. 

It is important to note that, in general, the end loads N x and N * , the inplane 
shearing loads N# and N#, the out-of-plane shear loads <?,' and Q* , the bending 
moments M~ and M * , and the shearing moments and can be known functions 
of 6. Also, the inplane shearing loads N& and N*, the circumferential loads N e ‘ and 
N e + , the out-of-plane shear loads <? e and Qi * die bending moments Af e ' and , and 
the shearing moments and can be known functions of x. The net distributed 
normal load q = q* - q~ can be a function of x and 6. 
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D. Application of the Method of Total Potential Energy to Cylindrical Panels 
Substituting the expressions for the strains given in eq. (12) into the energy 
expression for the cylindrical panel, eq. (5), and including the nine loading terms being 
considered, results in a rather lengthy but complete expression for the total potential 
energy, namely, 
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n(«»°) = 

\ f'i fi lit 


(',e-4)(Yrt^)]«d0<ir* 


+ / 9 *.V* (e) “°('T e ) Aie -/..!|^ (8) “1 t i- e h 9 

p P 

* /,' ! |W;(8)p;(-|,e)iM0 -/;>- (8)8 'H' e h 9 

p 

* / e *. 2 i N '» (e)v '('l’ 9 ) jMe 

2 2 

+ tit -/,'.!|< (e )^ + f' e r e 

♦ /;W|4^ -/.*i <? *' <e)w t*T e h e 

* ij.i w, wv -/„?| ^ 

L I 

L I 


+i 


*1 


- /J* fJ t qw°{x,e)MQdx 


(40) 


In the eq. (40), use has been made of eq. (11) and the fact that 


25 


dr = dz . 


( 41 ) 


and the integration with respect to r in the volume integral has been replaced with 
integration with respect to z. Grouping terms under single integrals and integrating with 
respect to z leads to the expression 


n(u °,v °,w °) = 

+1 


iUiL'AwY. 

*12 

- K-m/K 




♦I f 

* L-i 

2 l 

- qw°]RdBdx 

+ 

M;(e)p^|,ej - <(e)pf + |,e| 

+ 


+ 


+ 

(?,(6)w"|-|,ej - g;(e)w*|f|,ej 

* a { 

2 l 

- A£M»"fx,*.|jj 

+ 

M«*x)p£,-|) - ACWp£,*|) 

+ 

Ar,-(x)v‘|x,-|j - W,-(x)vjx,»|j] 

+ 

M.»P^x,-|j - «x)prfx,*|j 

+ 

C«(x)w*|x > -|j - Q,*(x)w*fx,+-|j J 


}RdQ 


( 42 ) 
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The stress resultants in eq. (42) are defined as 


H 

AT = f 2 b ojdz = i4 n e*+i4 12 eJ+yl 16 Yrt+B n ^+B 12 iCe+B 16 i4 -N? 

Z ’ 2 

N e = f \ a B dz = A l2 tl+A n el+A^y^B l2 i^+B 21 ^+B 26 Y^ i -N e 
Z ~~ 2 
+ H 

A^e = f _ 2 ff X x6 < ^ z = i4 l6 C Jt + ^26 E e + ^66Y*e + ^16 1C » + ^26 K e + ^66 K Jcfl ~^rt 

*- OOP 

M = f 2 h za x dz = B^isB^l+Brfrt+D^+D^+D^Krt - M x 

Jz—- 

B 

M e = f\ zojdz = fi 12 e%B 22 e“ + B 26 Y^e + ^ 12 < +I) 22^ +Z) 26 1< 4 ' M ! 

2 

*- o o P 

■ £ 2 W = B 16 e“+5 26 e2+B M Y*9 + ^i6 1 ^ + ^26 K e + ^« K *« * 


In the above 


B B 

' = / 2 * <&fc = / 2 „ (<? n e? + Q 12 el + 

J*- 2 I 

jj jy 

- f 2 h °e* ■ / 2 ff (^I2®r + ^22 c e + 

Jl ~T j*-_ 


2 

J5T 


J7 H 

/ 2 B X *& = f 2 B (Ql6*x + ^26 e e + 

J t .-_ J Z .- T 

*- - 

/__£ zoftfe = 2 * (Q u e' + <? 12 e e H 

*> 'l 


Qu1*)& 
Q^) & 

Qeelrt)^ 

h Qi 6 y F rt)^ 


M e F = [ 2 b za F dz = f 2 g (Q l2 t' + Q n 4 + Quirt) 

Jz—- J * m -j 

*2- *2 

M **{ m \ zx rt^ ■ j ]b (Qu z x + Q* e e + Q<*1rt) zdz • 


The superscript "P" denotes a preloading effect, and therefore, these expressions are the 
so-called equivalent preloading stress resultants. If the preloading effects are thermally 


induced, g g 

N x H ° x dz = f t __H {Qn a x + Qn a e + Qi6 a xe)^ Tdz 

2 2 

*- *- 

N * = N ! m /. 2 w ■ / 2 h Tdz 

2 *~~2 

N * " N * * L \ ^ = / 2 h (<?16«x + ^26 Ct e + ^66 a J A7 ’* 

2 * 2 

M x = M/ h | * Z o T x dz = f 2 h (Qu a x + Qu a e + Qi6 a xe) AT 
2 *“ 2 

- - 

M 9 = * f m \ Za J* “ / £ (<?12 a j C + ^22 a e + ^26 a J A7 ' 

*" 2 2 

^*9 ~ E f W J fj (fti6 a x + Q26 a Q + Q66 a X !&} A T zdz . 

2 2 

These expressions are the so-called equivalent thermal stress resultants. 

The notation 


(45) 


n = n («»*) (46) 

is being used in eq. (40) and (42) to emphasize the fact that the total potential energy is 
a function of the displacements, the superscript zero indicating that the displacements of 
interest are the displacements of the cylindrical reference surface. The governing 
equilibrium equations and consistent boundary conditions will be derived by examining 
the variations, or increment, in the total potential energy due to variations, or increments, 
in these displacements. Then, considering increments in the displacements of the form 
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(47) 


u°+eu{ ; v 0 +eV|° ; w°+zw° , 

the increment of the total potential energy will be of the form 

where e is a small parameter and the quantities u ° , Vj , and Wj satisfy all the 
kinematic requirements of the problem. 

The incremented total potential energy can be expanded using eq. (42) as follows: 
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n+An = { L-i /,.- ! ||K +AW .- w *')( e ° iAe ») * K*aw,-w,7«;>a.;) 

* Ary * (M,*AM r M^AK:) 

* (M,+AM,-M t j4+A4) ♦ (Afrf+AM^-M^icVAry 

- gjw °( jc, 0) + e w" (x,0)| iW0<fc 

♦ Ci [ w ;(4f‘l- e )* e “>t'l’ 9 )} ' ^4H e )-H e )l 

* M;(9){p^-|,ej«p;|-|,e| - <(e){p^|,e)-ep: ] [»|,ej ' 

* W rf (6) v*|-|,ej»ev 1 ^-|,6jJ - N>)jv^|,e)*tv^|,ej. 

* W b (x)|p^-|j + eP^x,-|| - ACf^^f )*<(*■ »f)[ 

+ W e (x)^x,-|jttv,"|x,-|| - W,*(x)|.*|x,t|j.ev 1 ‘|x, + || 

* M,-(x)|p^x,-|j + ep;|x.-|| - M,*(x^x,.|j +e p;|x,.|| 

* Oe(x>H'‘|x > -|Jtt>v 1 '|x,-|jJ - <? e , (x)|»*|x,»|j*ew 1 '|x.tl| dx . 

(49) 

It should be noted that the equivalent preloading stress resultants do not have increments 
because they depend only on the material properties, temperature, and initial 
displacements, not the displacements due to the applied forces. Subtracting eq. (42) from 
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eq. (49) leads to an expression for the increment in the total potential energy, namely, 


AE = 1 f} k / e A { K-iV/)Ae: ♦ AN x t° x * AN x Ac° x + 

+ A N e ee + AN e Aee + + Ai^Ay^ 

+ (M x -M/)Ak£ + AMX + AM x Ai^ + (M 8 -M/)Ak^ + AM^ 
+ AM 6 Ak£ + (m^-A^Ak* + AM^ + AM^Ak^ 
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The various strain increments in eq. (50) arise from the basic definitions of the strains 
and curvatures, i.e., eq. (13), by substituting the increments in the displacements, 


eq. (47), for the displacements in the original definitions. Such a 

substitution for e 

leads to 

. . d(u°+eu .") i „ , 

' P *.^ • 

( 51 ) 

Discarding the superscript zero for convenience and expanding, 


e +Ae =— +e— +-pf+ep B +-e 2 p 2 . 

x 1 ax ax 2 XXt 2 ** 

( 52 ) 

Substituting eq. (13) for t x , . . 

Ae =e — Up p )+— e 2 P 2 
* dx PlPx ‘j 2 Px ‘ 

( 53 ) 

which can be written as , 

Ae =ee +e^e 

X JCj Xj 9 

( 54 ) 

where e_ and e T are defined as 

*i *2 du. 

e = — Up p 
etc xPje ‘ 

( 55 ) 

e =— p 2 . 

*i 2 1 

( 56 ) 

In a similar fashion . 

ofv+ev.) w+ew, i 

JtdO * R + 2 <P#+eP ® |)2 

( 57 ) 

+t-^t-ip 2 .eBB +i-e 2 p. 2 , 
Rdd RdQ R R 2 P P#Pa » 2 Pa » 

( 58 ) 

or A 2 

Ae e =ee e 1 +C 8 e, ♦ 

( 59 ) 


where 
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(60) 


dv, w. 

Ce '~[^ae + T +p ® Pfli J 


E e =— P 0 2 

6j 2 r0 l 


(61) 


Likewise, the increment in is given by 




(62) 


dv ^1 du ^1 .no 

Y j» + Ay^ = + £ + + € + P-Pfl 

Yxe Yxe ax ax *a0 *ae 

(63) 




which can be written as . .7 

A YxB -8 »*e, 8 'rtj » 

(64) 

where 

y i^l + ^L + p p fl+ M ) 
Yjte > atx *ae PxPe > P6Px \l 

(65) 


Y x e 2 = Pjc 1 Pa 1 • 

(66) 

Finally, 

Ax,=eK Xi ; AK 8 =eK 9i ; Ak^-ck^ , 

(67) 

where 

#» 

I 

II 

* K ~ 

(68) 


a*.. 

Rde 9 

(69) 


K 

*®> ( ax jeaoJ 

(70) 
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( 71 ) 


In the above, use has been made of the definitions 

dw. dw. 

Px, ~ ~dx ’ Pe «" ~Rdd ’ 

The increments in the stress resultants, in terms of the strain increments just defined, 
are given by: 

AAfj, =i 4 n Ae x + ; 4 12 Ae e + ^i 6 Ay xfl + £ n Aie x + B l2 A Kq + B^Ak^ 

AN 0 = / 1 12 A e x + A22 A e e + A^Ay^ + B l2 Ak x + B 22 Ak 9 + B^Ax^ 

= ^16^ e x + ^26^ E e + -^66^ Y ^ + 5 16 Ak x + H 26 Al^ + AKjjJ ^ 2 ) 

AM X = B u Ae x + B 12 Ae e + B l6 A Yjrt + i>„Aic x + + D^Ak^ 

AAf e = B 12 Ae x + B 22 At e + B^Ay^ + D 12 Ax x + D^AKq + D^Ax^ 

= 5 16 A e, + ^26 A e e + *66 A Y* + ^> 16 Ak x + D^AiCg + Z^Aic^ . 

In order to separate the first, second, etc. variations of the total potential energy at a 
later point, it is convenient to expand the increments in the stress resultants in terms of 
the displacements and redefine those increments in terms of powers of e . Incorporating 
the definitions for the strain increments, eqs. (51)-(70), the stress resultant increments 
from eq. (72) are 



AAr x =^ n (ee X]+ e 2 eJ^ 12 (ee ei+ e 2 e aj ) + il 1< (eY xei+ e 2 Y xflj ) 

+ ^ii e \ + ^i2 eic e, + ^i«***e l • 

(73a) 

These terms can be redefined to give 



AN=tN +e 2 N , 

* z l *2 

(73b) 

where 

^Xj = An e Xl + A 12 e 8 l + A 16 Y x e, +B n\ +B u K * i + *i6 K *e, » 

(73c) 

and 

^x l = An e ^ + Ai2 e 8 2 + A|6Y*8 1 • 

(73d) 
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Expanding for future reference, 




.+— + 


Rde r 


du, w i 

la* P,P *>J 12 l 


e^o, 


+.4 


‘^a* *ae 


+ — - + B +B fl p r 

VM K * K ®I P ® K *l 


and 


ap„ . ap., fap ap ' 
^-aT^SS^-aT*™, 

^'(^„P,, J 4 ^ 12 P 8 l 2 ^, 6 P„P 8 ,) 


Using this procedure for the remaining stress resultants, 

+e 2 e J[j )+i4 22 (ee ei +e 2 e aj )+i4 26 (eY J( e | + e 2 Yj [ e I ) 

+B 12 eK x, +B 22 eK 9 1 +B 26 CK *e l • 


Redefining, 


with 


A N,=eW,« l N„ . 


w «,*Aij' Il + A B t,+A M Y ]iei tB, 2 K, i +S 2 2 Ke 1 t S M ': 1 ,o 1 


and 


^e J _ ^i 2 c x I + ^ 22 e e 2 + ^ 26 Y x ej 


Expanding for future reference, 


ap„ . aP 6 , „ ap,,) 

+Bl2 dx + “Rae^ 26 ! a* + *aej 


(73e) 

(730 

(74a) 

(74b) 

(74c) 

(74d) 


(74e) 
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and 


(740 

Likewise, 

A ^e = ^i6( e S +c2 S )+i4 2«( ce e 1 +c2c e 2 ) + ^66( c Y J( e | + e 2 Y Jiej ) 

(75a) 


+ ^ 16 elC *, + ^ 26 eK flj + ^ 66 CK *e l * 

or 

A ^e=eA^ i+ e 2 A^ , 

(75b) 

where 

16 C x, +i *26 C 8 I + -^66Y ifl, + ^ 16 IC *, + ® 26 K e | + ^66 K *e, 

(75c) 

and 

A ^e J = ^i6 e * l + ^26 8 e J + ‘^«6Yie J • 

(75d) 

Expanding, 




dx u Rd6 dx Rddj 

(75e) 

and 

w -,=(^«p, 1 2 4^p. 1 2 ^«p,,p. 1 ) • 

(750 

The increments in the moments can be similarly defined, namely, 



AM x =B u (t t Xi *t\) +B n (t e 6i +e 2 Yjej ) 

(76a) 


+D n eK, i +D 12 eic ei +i> w eK <6i , 

or 


(76b) 

with 

= ^n e X[ + -®i 2 c e 1 + ^i6Yxe 1 + ^ > ii K x 1 + ^i 2 1c e 1 + ^i« K *a 1 

(76c) 
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and 


^x, _ ^n c j(j + ^i2 e e 2 + ^i6Yj(e J » 


where 


and 


M x t =B i 


+B 


1<H 


(du. 

— -+P P 

dx PxPx ' 

^♦±U Mt+W ,' 

dx RdQ * e ‘ 6 1 


( dv, w. S 

+fil2 [^0 + ~^ +PePe '> 




(dK 3 P 


'll 


dx 


'12 


Rde 


i« 


i, i 


^ dx Rdd , 


2 i ^*i 2 + 2 ^ fl » + ^ 16 ^ Jt i^ e i] 


In a similar manner, 


A Af e =B l2 (e +e 2 e ^ +B Il (e e fli + c 2 e 6j ) +B 26 (.e +c2 Y J[ e 1 ) 

+Z) 12 eK * l +I) 22 eK a 1 +I) 26 elC i fl l » 


or 


with 


and 


where 


and 


AM e =eM 6i +e 2 M 6i , 
^ 6 2 = ^ 12 8 ^ + ^ 22 e 8 1 + ^ 26 Yx 8 I » 



(76d) 

(76e) 

(76f) 

(77a) 

(77b) 

(77c) 

(77d) 

(77e) 
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Finally, 


where 


with 


and 


^-(|fi,^|» 2 2P., 2 *»*P, 1 P. 1 ) • 

AAf *e =fi i6( 6 V +t \) +B u(c\ +c2c e J ) +5 «6( c Y Jf e I ^Y^,) 
+D u>*\ +D 26 tK e [ + D 66 eK x e l » 

^ = ^16 e x, + '®26 e 0 1 + ^«Y J a 1 + ^16 K * i + ^ 26 K 9, + ^66 K xe, 




Expanding for future reference 


^*e. “Ad 


^“i ) f dv, w, > 

— ±+PA +fi J— L+ - i + P«P« 
a* x x> ) 24 \/j80 * 6 e «J 

uaj^u**) 

> *\ ax *aej 


ap„ 5 Pe, 
16 ax M *ae 


and 


M ^=(|*,eP,, J 4^P 8 , 2 *»«P Il Pe 1 ) 


(770 

(78a) 

(78b) 

(78c) 

(78d) 


(78e) 


(780 


With the strain and stress resultant increments defined and expanded, the definitions can 
be substituted into eq. (SO). This results in 
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+ 

+ 


ad - \ 0 f e '.i (K-<K** 2 s) * («v e X)*. 

Z 2 2 

( e V e2N J( c v £ \) + 

(tN^N^t, * (tiV e n^«,n l e,J 

K-^X c Y A ^ 2 Yje2 ) + (eJV'Xe,)^ 
(eJV^+e a N^eY A +«%J + K _A 0 £K *. 

{tM^t 2 M h ) k x + (eM^+e^JcK^ + (AVA//^ 
(cA^+e’AfJi c, + (eM^e'M.Je^ + (M^-M^tK^ 
(eAf^+e’Af^K* + 

■{0{h><4 e )- w;<9) “'H e )] 

♦ m;< 6) p;0.e) - p'0 ,e ) 

- Aimvd-^j - N»v,^|.ej] 

♦ M*(e) p;/-|,e) - *C(6) 0f,e 

♦ <?;(0)w,“f-|, e j - <?,*(« 0|.ej 

Lh {fw44) - ^>4 ,+ i)l 

♦ M*(x) p^x,-|j - M^x) P^+f )] 

+ N e "(x)v 1 °|x,-|j - <(x)v^x, + || 

+ M e (x) - W «W P e( X,+ l)] 

♦ <?e‘C*)w l *fx,-|j - <?;W w i°( x ’ + f )] } dx 
Expanding and regrouping in powers of e leads to 


(79) 
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An = 


+ 


\ Ui Ci 4 n .-n.% * v * MX * v* 

46 2 2 L 

* - W„,Y* * (*f,-*X * M^ic, 

* fc-*X * M„k, * - «*,*} 

* « J {fc-«X * *v* * w -,s * K-*X 

+ ^ej e e + ^o^e, + (^d _ ^U,)Y*e + + ^*e,Yie, 

+ A^k, + + AVe + ^e,^, + Af^ + A/***] 

+ e3 Ks + N **\ + V®, + V®, + N *?a* 

+ + + Af^J 

* |t?| {[ w /( 0 )“i(-|. e ] - N;(0)«^|,e)] 

* «;<«) p;|-|. e ) - «;(8 ) p;/*|, 8| 

* Wv,‘(-|.ej - Aa»)v,*(*|,e| 

* *Q0) p;(-|,e) - m» Pjf|,e| 

* <?;te)H' 1 j-|,e| - <?;(8)H.,^|,0j] 1 

e {A { A' ) ■ ^“^'"f)] 

* Mi w - w^x) 

* JV."(x)v,‘jx,-|j - W«(x)v,‘|x*'|j] 

* m,(x) p;|x,-|j - Mia p;|x,*|| 

+ C«(x)w,‘|x,-|j - <?,’(x)w l 'fi,->-|j] •< 


RdQ 


( 80 ) 
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The increment in the total potential energy can be written as 


An=en i +e 2 n 2 +e 3 n 3 +e 4 n 4 


( 81 ) 


The quantities ell^ e 2 ^, e 3 ILj, and e 4 H 4 are defined to be the first, second, third, and 
fourth variations of total potential energy, respectively. The equilibrium conditions for 
a cylindrical panel are obtained from the condition 

enj(Uj,Vj,Wj)=0 - n^.v^-o , (82) 


where the notation indicates D, is to be stationary with respect to the displacements 
n Vj, Wj . These displacements are the variations in the equilibrium displacements. 
The second variation is used to examine stability of the equilibrium displacements. 
According to the Trefftz stability criterion, transition from a stable equilibrium 
configuration to an unstable one is characterized by 

Se^KpVpWj-O - AH^iipVpWj-O . (83) 

This states that the second variation of the total potential energy should be stationary with 
respect to variations in u v v,, w 1 . 

Interest in the present study focusses on equilibrium rather than stability. Hence, the 
first variation, , shall be studied further. The higher order variations, however, will 
not be discussed beyond this point. 

The first variation can be identified with n, and that quantity is given by 
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n rjUi L-i l K rK%^(N,-N,y *n t, 

^22 

*( M e ~M^ I * M e, K »*( M «i “ M ieK», * M a, K rt- < l w ° ) *** 


fill jp«)“.'(-|.e) - <(e>“.*( + f.e)J 

* K( 0 ) p;(-|.s) - «») .e)j 

* ^(e)v,“|-|,ej - w»v,"|»|,0j] 

* w«(6) p 5(-|- 6 ) - p «(*f’ e )] 

* <?,-(«) - <?;(e)^|^]Uje 

f,.!i { s m “‘”(^'"1 ) " '*».« )] 

mm p:,(*.-f) - mm p^*|) 

Wwv,*|x,-|j - W,*(x)» 1 ‘|x, + |j] 

mm (>;,(*,-{) - MM ps^^I)] 

<?e(x)>*i'(*.-|] - Oe(x)w,‘|x,*-|jj 


(84) 


A more useful form of the first variation can be obtained by substituting for N , W e , 
^*6j » * M o t » 311(1 M to l from eqs. (73e), (74e), (75e), (76e), (77e), and (78e). If this 

is done and the various terms in this expanded form of IIj are regrouped, the result is 
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.i .1 

n, = L i /.- 2 i 

2 2 

* ) BdMx 


♦1 [\ 
4 1 

w,'C0) “»”( -f > e ) ’ « e >“>"(*|- 6 )] 

+ 

m;< 9) p;(4 e ] - « e > p“,( + 2 ,e ) 

+ 

w^e)v 1 'j-|. e j - N»v 1 "|4 6 jJ 

+ 

M rt (0) p;,(-| ,e ) ' M>) p; H’ e )l 

+ 

<?;(9)4'i’ e ) " <? ' (e)w, ”(*2' e )] 

- 4 ( 


+ 

mux) p°41)J 

+ 

- W.*(x)v l *|x, + |)] 

+ 

Mh(x) " M 9 (x) P e i[ I,+ "2^j] 

+ 



Note the quantities A//,... Af* have disappeared, as has the factor of 1/2 in front of the 
two-dimensional integral. Substituting the strain-displacement and curvature-displacement 
relations, eqs. (55), (60), (65), and (68)-(70), into the above, the first variation takes the 

form 
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This is one of the fundamental forms of the first variation for determining the response 
of a cylindrical panel. This form can be used directly in approximate schemes such as 
the Rayleigh-Ritz method. Further steps, however, are required to obtain the equilibrium 
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equations and associated boundary conditions which are of interest here. 


E. Application of Integration By Parts 

To determine the equilibrium equations and the associated boundary conditions, 
differentiation of u lt Vj , and w x with respect to the spatial variables x and 6 must be 
eliminated. This is done by applying integration by parts to the various terms in 
eq. (86). This procedure follows, the results being given on a term-by-term basis: 
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fifth term 


.2 


♦2 


f \ f %W e P e P 0 M 6 dx = f \ ( 2 . N 9 — ^ RdQdx 
J x ,-± Je.-l er » r0 i J*=-2 J e «-2 6 jma bar 


RddRdd 


sixth term 


seventh term 


■ A * 


2 

■ L *2 


- f 2 L f 2 p — K— Vi MQdx 

J *~f Je-I *Rdej 1 

1,1 /i "-It*®* ■ /.i pv-Cl *• 


+2 


r 2 r*! 

-Ilf B — — v.iWedlx 
J,-| Je— 2 ar 1 


♦2 


a /A ■ a mCI * 


♦2 


♦2 


- a a 


3AL 


'*~f J «-£ RdO 1 


u.RdQdx 


( 90 ) 


(91) 


(92) 


eighth term 



♦2 


■ f' f> Hjgu MBA 


a* jwo 




( 93 ) 
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ninth term 



N^£,RMdx - [\ /j. 


'Rdd dx 

L 


■ L'{ (»-s-R " 

- ri lit IKSi-'”* 


( 94 ) 


tenth term 



(95) 


The last term on the right can be expanded further by substituting the definition for , 
from eq. (71), into eq. (95), namely, 



♦1 


RdQdx = f 2 l L 2 i 

fa x, Jjc*-— J 0-— C- 


dM x dw 1 
dx dx 


RdQdx . 


(96) 


Using integration by parts yields 



The tenth term can, therefore, be written as 
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eleventh term 


♦2 


ap 


♦ k 


.♦2 


LU LU - IJk ( M MJ * 


♦2 


(99) 


aM n 


- / f % — -p e /eded^c 

j e -2 K «, 


This can be expanded further by substituting the definition for P 0j , from eq. (71), into 
eq. (99), namely, 


f‘i bm, 

*ae 


\RMdx = £ 


J?36l?a0 




( 100 ) 


Using integration by parts and substituting into eq. (99), the eleventh term becomes 


♦2 


C2 t*2 d K 
life Af e — —RdOdx 
J x .-± Jo-2 « DM 


♦2 


= - £_l (w. 
*2 / 

- - 


dw. 


l 


/?ae 

♦2 


-2 

2 dx 
-2 

2 


A A 


Rd6 RdQ 


(101) 


f 2 f 2 

- / 2 f 2 S-w./MBdr . 

J,-| Jo-2 *2^ i 
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twelfth term 


♦i 

2 


Ui L\ M » 


ap, 


dx 


-RdQdx 



- s2k /. 


2 

-4 


dx 


P 0 RdQdx 


Using the definition of p 0j , in eq. (102), 



Expanding the second term on the right using integration by parts, 



dp e 

— -RdQdx = - 
dx 







r 2 f'i £a* w j UUt ' 

Jx~ | Je-{ RdQdx 1 


Applying integration by parts to the first term on the right, 


(102) 


(103) 


(104) 
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(105) 



Collecting terms evaluated at -L/2 and +L/2, the first term on the right in eq. (104) 


becomes 



(106) 


The twelfth term can thus be written 



(107) 
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Expanding the first two terms on the right for future use, the twelfth term becomes 


*k *1 ap n 

f \ ( \ M^—MQdx 

J x ..i Je~i * dx 


.1 r->w *- \ am lb-+i 


/. 


2 


'?*»♦/ \ 
L J jr— — 




+ I . I — — w. 

-f [*ae 1 

*k *i. &M 

_f 2 f 2 mdx 

J*~4 Je-4 j?aea* 1 


ebe 


NLl* 


(108) 


Thirteenth term 


+- +-& 0B f*k P“*4 

Ui Lh - Lt p-mU* 


- Li /, 


+1 

2 


dM, 


J«— | RdQ x 1 


2 -V x RdBdx 


(109) 


Substituting the definition of p Xj into the first term on the right above and integrating by 


parts: 


.k 




.k 

2 


dw 


"-irlt’-j* 


- /:!t[4*f)5(‘*l)]* • 4hH)3K] 


dx 


(110) 
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Substituting the definition of P Zj into the second term on the right of eq. (109), 

. dw. 


r 


2 

r»-— 

2 

♦1 


♦1 

f 2 . 

8M. 
19 1 

2 

Rdd 


Y* 

Rdd 

■i- 




( 112 ) 


Substituting the results of eq. (Ill) and eq. (112) into eq. (109), the thirteenth term can 
be written as 



(113) 


Expanding the first two terms on the right for future use, 
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♦I 


3P, 


L± /..!± * 

2 2 


1 /eae 


- m *B >+ f H*f ,+ 2 ) + 

Ht3- 


/ \ 


4 [aw: 


*^ + r*! I« 

1 Je— 


4 hae 


♦i 


- Ut /. 


4 ato 


2 e~j Jiaecbc 


— y v.Rdddx . 
1 


(114) 


The terms involving Af*. and in the line integrals in eq. (86) must 
be expanded in terms of Wj and integrated by parts also. These two terms follow. 


Term involving and M# in line integral from g=-fl/2 K> 8=±8I2l 
Substituting the definition of P 0j into this term and integrating by parts, 

hit [ M -( 0 >p.,(-|’ e ) ■ M ^ e)0 »,(*|- e | Rde 




(115) 


Expanding the first and third terms on the right-hand side, 
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(116) 



Term jpvplying and in line integral from x=-L/2 to x = +L/2: 
Substituting the definition of into this term and integrating by parts, 

Vi [*w., (*.-{) - 

■ ' Vi * Vi ) A 

■ ■ * A f ^ w 'Kh 


A/^Wh', x,+-|j 


*-4 


m ~2 


3 A dU^x) 

L h-i 


dx 


f*!)* ■ 


Expanding the Erst and third terms on the right-hand side, 


(117) 


54 



(118) 



Substituting eqs. (87), (88), (89), (90), (91), (92), (93), (94), (98), (101), (108), 
(114), (116), and (118) into the expression for the first variation, eq. (86), and 
combining the boundary terms (i.e., integrals with respect to 6 and integrals with respect 
to x), results in 
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For the first variation to be zero, each term which is a product involving an increment 
in the displacement, in each of the integrals, must be zero. The Euler equations come 
from the two-dimensional integral and the boundary conditions come from the line 
integrals along the x and 0 edges. Thus, the three governing equilibrium equations are 



(120a) 


(120b) 


(120c) 


Using the first two equations in the third one, the three equations can be written as 



£&♦££- o 

(121a) 


dx RdQ 



w* + =0 

(121b) 


dx RdQ 


#M X 

"to* 

„ A x , Sw 

* 2 «aea i rW ‘ex 2 

(121c) 


Rdddx R 2 dQ 2 R 



The variationally consistent boundary conditions at the x L/2 and x +L/2 edges of 
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the cylindrical panel are 
along the x = -L/2 edge: 


i) 

N x = 

n; , 




or 

u must be specified. 

ii) 

= 

A(*e * 




or 

v must be specified. 

iii) 

N& 

+ AT* 

dw 

dM x 

+ 

8M„ 

+ 2 


- dM * 

Q x + * , 

X dx 

Xu 

/?ae 

dr 

«ae 


Rdd 







or 

w must be specified. 

iv) 

M x = 

w; , 




or 

dw 

— must be specified. 
dx 

2LE. 

+L/2 

edge; 






i) 

A^t = 

at; , 




or 

u must be specified. 

ii) 

A^e - 

A&, 




or 

v must be specified. 

iii) 


+ ATa- 

dw 

dM x 

+ 

+ 2 

= ( 

* dMl 
?; + * , 

dx 

XU 

f?ae 

dr 

/we 


* RdQ 







or 

w must be specified. 

iv) 

M x = 

a/; , 




or 

dw 

— must be specified, 
dr 


(122a) 


(122b) 


The variationally consistent boundary conditions at the 0=-/3/2 and 8 = +/3/2 edges are: 
along the $ = -812 edge; 


i) 

ii) 

iii) 


iv) 


AT* * N* , 
- A/ e . 


or u must be specified, 
or v must be specified. 


N a 


dw 
'Rd e 


dw SMI 


e 


dr Rd 0 


>r vrr 

AT.— + 


+ 2 - 


8AT 


= <?•' 


5Ma 


M e = Af e " , 


& ' w dr 

or w must be specified, 
dw 

or — must be specified, 
dr 


(123a) 
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along the 0 = + 13/2 edge: 


or u must be specified, 
or v must be specified, 
dMi> r 


i) ** = K , 

ii) > 

„ 8w , „ . 3M 0 (123b) 

6 l?ae *dx Rd6 dx dx 

or w must be specified. 


iii) •‘’e 


iv) M e = M e + , 


or must be specified. 

dx 


Note that the expression for the first variation of the total potential energy also involves 
four non-integral terms. These are "comer" conditions. Therefore, 
at the corner (x = +L/2. & = ± d/2) ; 


2M a = , or w must be specified , 

at the corner (x = +L/2. 0 = -8/2) ; 


(124a) 


2 , 


or w must be specified , 


(124b) 


at the corner (x = -L /2. fl = +fl/2) : 

2^ = M* + ML , or w most specified , 
at the corner (x = -L/2. fl - -fl/2) ; 


(124c) 


(124d) 


2M rf = + A/a, , or w must be specified . 

At this point, the problem can be specialized to that of a complete cylinder, that is 


P =2 71 -*■ 


-A = -« +— = +n 

2 2 


(125) 


Because the cylinder is complete, the response is a continuous function along 0=x, and 
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the displacements at the -j8/2 and +/S/2 "edges" are unique and therefore can be 
interpreted as being specified. Equations (123) do not have to be explicitly enforced. 
This also holds for the "comer” conditions. However, a by-product of the eqs. (123) is 
an expression for the out-of-plane stress resultant, Q 0 , which could be used as a check 
of a calculation of the interlaminar stress component x 6r . Specifically, 


<?e 


N a 


dw 
Rd e 


+ AT. 


x6 


dw 

dx 


+ dM, 

RdQ 


dM « 

+ — — . 

ax 


( 126 a) 


This shear stress resultant compliments the other out-of-plane shear stress resultant, Q x , 
which is given by 


Q, 


- N x — *N„— 
1 dx 10 RdQ 


dM dM~ 

+ + — . 

dx RdQ 


( 126 b) 


This out-of-plane stress resultant will be used as a check of a calculation of the 
interlaminar stress component x xr in chapter VI. 

Attention now turns to specializing the equations for the condition of axisymmetry. 
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m. SIMPLIFICATION OF THE EQUILIBRIUM EQUATIONS 
DUE TO THE CONDITION OF AXISYMMETRY 


The derivations which follow focus on complete cylinders (i.e., (3=2r) subjected to 
an axisymmetric loading and responding in an axisymmetric manner. In particular, the 
loading considered will be applied axisymmetric end loads N S (Q)=N X and N X (Q)=N X . 
Therefore, the relations of the previous chapter will be simplified by 

iQ=Oand-^i=^, < 127 > 

36 dx dx 

( ) being any response quantity. Under these conditions the kinematic relations simplify 
considerably. Specifically, eqs. (9) and (13) become 


p;=~; p:-o 


dx 


du 


1 no* 




dv* 


8 x = — e e p - ’ Y,e & 


dx 


(128) 


d P 


< = kS=0;4,=0. 

dx 


The equilibrium equations, eq. (121) simplify to: 



dx 


0 


dx 


= 0 


(129a) 


(129b) 
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(129c) 


d 2 M x 

dx 2 


+ N. 


dx 2 



The accompanying boundary conditions at the x=-L/2 and x = +L/2 boundaries are 


at x = - L/2 


i) 

N x = 


or 

ii) 

Ifl 


or 


dM r 

dw 


iii) 

X 

dx 

■ q - ■ 

or 

iv) 

ii 

af 

u; , 

or 


u must be specified, 
v must be specified, 

w must be specified, 

— must be specified. 
dx 


(130a) 


at x = + L/2 


i) 

= K . 

or u must be specified. 


ii) 


or v must be specified. 


iii) 

dM „ Aw . 

-3T - q: • 

or w must be specified. 

(130b) 

iv) 

II 

dw 

or — must be specified. 

dx 



Note that for a complete cylinder, the boundary conditions at B=-\ 3/2 and 6= +/3/2 are 
meaningless. Equations (129) and (130) will be the focus of the remainder of the this 
chapter and the following chapter. In the next section, the solution of these equations 
for the case of a known axial end load will be derived. 
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A. Solution of Equations for the Case of Axial End Load 
The first equilibrium equation, eq. (129a), integrates to 

N = constant. 

jr 


(131) 


Since the axial load is known at the ends of the cylinder, this constant is the applied end 
load. It will be referred to simply as N. The second equation, eq. (129b), integrates to 

V- = another constant . ^2) 

JTO 


In this study there will be no external torsional load applied explicitly to the ends of 
the cylinder. Rather, the tangential displacement of the end will be specified. 
Enforcement of the tangential end displacement may well require a torsional load, S , on 
the end. This torsional load will be solved for as part of the analysis. Thus, as a result 

of eq. (132), 


N a = S (a constant). 


( 133 ) 


To solve the third equilibrium equation, eq. (129c), it is convenient to express all 
quantities in that equation in terms of the midplane displacement w°. Here, the stress 
resultants, eq. (43), with the preloading condition being thermal according to eq. (45), 
will be used. Only cylinders which are balanced, i.e., those with A 16 , A^, and N# 
equal to zero, will be considered. This is the situation of many laminates of practical 
interest. Such laminates are called balanced laminates. Expanding 

N x , N e , N a , and M x in terms of w°, t° x , and y*e results in 
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(134) 


N x = A n t° x + A u — - B u ^- - N x = N 

x n X 12 R 11 &2 


^9 “ ^ 12 C x + ^22 


o - w° „ dV' 


- B 


12 ' 


<fr 2 


- IV/ 


^ - -4«Y- - - S 

M . - + i, 2 -y * «.,Yi - - Ml . 

Solving the equation for N x for t° x , and the equation for N a for y# yields 


0 

E , = 


*11 


N + N x - A n — + B u =-^~ 
x 12 * 11 A 2 


(135a) 


and 


0 

Yrf = 


( *16 d 2 ** + 

A* 


, ^66 "66 / 

Substituting these results into the equations for and M x in eq. (134) yields 


(135b) 


"e = 


A 

-4, V 


L *»^ 

^22 " 

24, J 

/? 

h‘ 24,, J 


dhv* 


d x 2 


>♦*/)-< . 


*11 


(136a) 


and 


M x = 


* 12 - 


*,1^,2] 

w° _ 

f lj2 > 

n "16 "11 

d\f° 

24,, J 

R 

11 7" T” 

^ ^66 a 11 ) 

dx 2 




Ml 


66 


(136b) 


Substituting eq. (136a) and (136b) into the third equilibrium equation, eq. (129c), leads 
to the equation governing w° : 
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(137) 


li 


Bn 

4' 

d*w° 

^ii 


dx 4 


22 


41 


■lU 


_ 

/e 2 


* 


^ 11^12 />^12 
( 2 vr * 



jV 0 

<ix 2 



This is a linear differential equation with constant and known coefficients which 
depend on material properties, equivalent thermal loads, geometry, and a known applied 
axial load N. Hence, this structurally nonlinear problem results in a mathematically 
linear problem. Note, however, that the coefficients vary as the applied axial load N 
varies. This will be addressed in a later section. 


B. Specification of Boundary Conditions 

Attention now will be focused on the boundary conditions. A statement has been 
made regarding N x . This statement satisfies one of the four boundary conditions, 
namely, eqs. (130a) and (130b), i. For the problem to be properly posed, the three 
remaining boundary conditions, eqs. (130a) and (130b), ii, iii, and iv, must be satisfied. 
Within the context of the admissible conditions of thermal preloading, general balanced 
laminates, and no explicitly applied torsional loads, three physically plausible boundary 
conditions can be imposed on the ends of the cylinder, namely. 

1 - lubricated boundaries; 

2 - simply supported boundaries; and, 

3 - clamped boundaries. 

In the remainder of this work, it will be assumed that the same boundary conditions 
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will be enforced at both ends of the cylinder. Hence, the discussions which follow will 
focus on the x = +L/2 boundary only. Since the same boundary conditions will be 
enforced at either end, the radial response of the cylinder will be assumed to be 
symmetric about the x=0 plane. Therefore, the radial displacement, w °(x) , of the 
cylinder is an even function of x. 

For lubricated boundaries, the ends of the cylinder are free to rotate about the 

cylinder’s centerline, implying that = S = 0. Also, the out-of-plane shear force Q x 

and moment M x at the ends are zero. For simply supported boundaries, the tangential 

displacement v° and the radial displacement w° at the ends of the cylinder are specified, 

and the moment M x is zero. For clamped boundaries, the tangential displacement v°, 

dw° 

the radial displacement w®, and slope —j- at the ends of the cylinder are specified. 
The terminology ’lubricated’ comes from the fact that there is no restraint on the 
tangential or radial displacement or slope at the cylinder ends, as if the ends of the 
cylinder were being pushed together axially with perfectly lubricated plates. 

From eq. (130), the formal boundary conditions for lubricated boundaries are 

= 5 = 0 ’ (138a) 

2 

= 0 , (138b) 

2 2 
and 


dx 


dw c 
+ N— 

xm+i. dx 
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(138c) 



For simply supported boundaries the conditions are 

v° I l is specified , 

w° I l is specified , 



For clamped boundaries the conditions are 

v° | l is specified , 


(139a) 


(139b) 


(139c) 


(140a) 


w° l is specified , 
**♦ — 


(140b) 


and 



is specified . 


(140c) 


In order to conveniently impose these boundary conditions, they will be expressed in 
terms of the radial displacement w°(x) and its derivatives. Therefore, the tangential 
displacement v°(x) is required as a function of w e (x) and its derivatives. From eq. 
(128), it follows that 
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^£-ri - v"(i) = f(r*(x))dx * c v , (i4i) 

where C v is a constant of integration. Substituting eq. (135b) into eq. (141) and 
integrating results in, 


V%0 = ♦ J-X * c„ . 

4 * * ^66 


(142) 


For convenience , and since we expect v°(jc) to be odd since w°(x) is even in x, 


C v = 0 


(143) 


In the following, each of the three sets of boundary condition equations, eqs. (138), 
(139), and (140), will be presented in terms of the radial displacement w°(x). 

L Lubricated Boundary Conditions 

Despite the desire to write the boundary conditions in terms of the displacements, the 
first of the equations, eqs. (138), which specify lubricated boundary conditions, i.e., eq. 
(138a), can be represented most simply as 

5 = 0. (144) 


The two remaining equations, eqs. (138b) and (138c), can be expanded to yield 


*16 
D — — 

r 2) 

d 3 w° 

+ 

^11^12 

V 

dw° 

11 A 
A 66 

~ A n\ 

dx 3 

2 



! ^ 


-1 


= 0 , (145a) 


and 
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(145b) 


( 


B 


2 \ 


D - _li 

11 A A 
A 66 A \\) 

B, 


d 2 w° 


dx 7 


H 


- -±(N+N x T ) * Af/ 

^11 


^ 11^12 _ ^ 12 ^ 

*, 

^5 = 0 . 
*66 


W °\ x ..k 


Equation (145b) can be simplified by enforcing the first equation, i.e., by setting S - 0 . 
Due to these simplifications, the case of lubricated boundary conditions results in the 

following two equations, 


ft 2 

n B 16 

B 2 ' 
g ll 

d 3 w° 

4 

^11*12 **12^ 

dw° 

° n A 
V A 66 

*ll y 

dx 3 

L 

X “*~2 

A n R R j 

dx 


-I 


. o , (146a) 


and 


*>n" 


*1 

D 2 ^ 
^11 

d 2 * 0 

4 

*^11*12 




~ A J 

dx 2 


1*11* 

R J x 2 

(146b) 

- 

2h»-< 

^11 

|<M,' 

= 0 . 




2. Simply Supported Bou ndary Conditions 

Expanding eqs. (139) in terms of w° yields 


L = 
2 


*16 dw° 


' x -*i Am dx 




is specified , 


(147a) 


w ° I l is specified , 
x ~*j 


(147b) 


and 
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(147c) 


*>n- 


66 


i>2 ' 
"11 

d^w 0 

^11 > 

dx 1 


2 


* J "• 


- — (W+W x r ) + M/ + —5 = 0 . 


‘u 
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^ Clamped Boundary Conditions 
Expanding eq. (140) in terms of w°(x) yields 


But dw* 
V ° £ = — — 

'■*2 ii * 


* -(+-) 

2j 


is specified , (148a) 


w°\ x _ l is specified , 


(148b) 


and 


dw ‘ 


dx 


2 


is specified 


(148c) 


C. Solution of the Governing Equation for w°(x) 

Equation (137) will now be solved for w°(x). The complete solution to eq. (137) 
consists of homogenous and particular parts, i.e., 

W ( (x) = W^(x) + w^fCx) . (149) 

The particular solution is simply the right hand side of eq. (137) divided by the 
coefficient of w°, or 
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(150) 






W. 


^ 11^22 ^12 


/wrt. 


where the notation w denotes the particular solution. Note that it is not a function 
of x. 

The homogeneous solution is of the form 


O) - 


(151) 


Substituting this assumed form into eq. (137) results in the characteristic equation: 


n ah ah ' 1 

u 'i,, 


l * j 


^22 ” 


n YL 

A nJ 


= 0 , 


or 


R* A ' 
n 4 _B 2 1611 

^ll-^n "it 

A 66 


X 4 * (2B n A n -2B a A n -A n RN)X 2 * 


•^ 11^22 ~ A 12 


= 0 . 

(152) 


There are four roots to this equation, namely 


*iam * 


t 

\ 


RN ♦ 2A lt B u - 2j4,20,|)±^ 

pA a B n -2A n B u -A u RNf-* 

B 1 A 

n a b 2 “ “ 

v u A li "II A 

i -*—l 

|(^1 1^22 

2 

2 

J 

R 


(153) 


Though it is not totally obvious, there is an interesting character to the roots given above. 
This is due to the dependence of the roots on the level of the applied axial load, N. The 
character of these roots can be examined by studying X 2 instead of X , i.e., 
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n, = 


(A u RN+2A n B u -2A a B u )±' 








(154) 


The first important character to observe is that the discriminant in eq. (154) will be zero 
for a certain level of applied axial load N. This load will be denoted by N * , given as 


AT = 


2 


A U R 


N 


(^11^22 ^12) 


\ A 66 ) 


+i4 12 B n A n B n 


(155) 


In general, eq. (155) will yield one positive and one negative value for N* . In this 
investigation, attention will be focused on compressive loading of the cylinder and, 
therefore, on the negative value of N * . Hence, N* will be given by 


N* = 


2 i 


A n R 


(<4ir^22 ^ 12 ) 


V A 66 ) 


+A l7 B n A n B l7 (. 


(156) 


For varying values of N, the roots (X 2 ) u have the following characters: 
1) For 0 i N > N * , the roots (X 2 ) u are complex conjugates given by 
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(^ 2 )u = ± 




<V„-».’.-^ * 


The roots X are thus of the form 


*UA4 * " ±a±, P 


(158) 


2) For N = N* , the roots (X 2 ) u are negative repeating real roots given by 


{A n RN*2A u B a -2A u B u ) 

<*>“ ' , ' 

— r^r 


(159) 


Therefore, the four roots X are two pure imaginary repeating roots of the form 


^u3,4 ■ ± - ** 


(160) 


3) For N < N* , the roots (X 2 ) u are negative distinct real roots given by 


(x*) w * 


(^n ®ia~^ia®n) ± ^ 


r >1 ' 

n ^ -b* - “ " 

i SLJ 


2 

i *iV*nl 

it 


(161) 


The roots X are thus of the form 
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= *V^U = ±l ‘Pl» ±*Pa * (162) 

Because of these three different forms for the roots of the characteristic equation 
(152), depending on the value of N relative to N* , the functional form of the x 
dependence of the homogeneous solution depends on the value of N. Therefore, the 


value of N defines the shape , as well as the amplitude of the deformed cylinder. For 

dhv 
dx : 

amplitude, not the shape, is dependent on the value of N. 


the linear problem where the N x term in eq. (129c) is not present, only the 


The functional form of the homogeneous solution, eq. (151), is as follows: 

For 0 zN>N* , from eq. (158) 

CW - + i 4 3 e ( -« + *>* + i4 4 e ( -*' ip)x . (163) 

For N = N* , from eq. (160) 


*W(*) ■ (*i + A 2 x)e l > x + (A 3 +A 4 x)t^ x . < 164 > 

For N < N* , from eq. (162) 

w^ ix) = A^'* + + A 3 c-*' x * A 4 <r** . (W5) 

Combining these homogeneous solutions with the particular solution, eq. (150), and 
considering only the portion of the solution that is symmetric about x=0, the three forms 
of the solution of w°(x) are: 

For 0zN>N* , 
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w°(x) = F cosh(ax)cos(Px) + G sinh(ax)sin(Px) + 


■^ 11^22 ”^12 


( 166 ) 


For N = AT, 


w°(x) = Fcos(px) + Gxsin(px) 


+ R j^n N 9~ A \2( N ~ N x)) ( 167 ) 


^ 11^22 ~^12 


For N < N*, 


w°(x) = Fcos(P x x) +Gcos(P 2 x) 


•^11^22 -i ^12 2 


( 168 ) 


The constants F and G can be determined from the application of the boundary 
conditions, eq. (146), (147), and (148), for lubricated, simply supported, and clamped 
boundary conditions, respectively. 


D. Solution of the Governing Equation for u°(x) 

The remaining displacement variable u °(x) can be obtained from the definition of the 
midplane axial strain t° x (x), given by eq. (135a), in terms of the material properties, 
cylinder geometry, and equivalent thermal loads, along with the solution for w°(x) and 
its derivatives. From eq. (128), it follows that 


du° (x) 1 f dw°(x) Y 

dx ‘ xU 2l dx 



(169) 
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where C u is a constant integration. By substituting eq. (135a) into eq. (169), 


du°(x) _ 1 
dx A 


u 


dx 7 


) 1 

\2 

dw°(x) 

j 2 

dx J 


( 170 ) 


This expression involves constants and the closed-form solution for w°(jc) and its first 
two derivatives. Therefore, this expression can be integrated analytically as 
characterized in eqs. (169), i.e., 


u°(x) = 


— (n + w/)L - f 

V '] [Ai*JJ 


w°(x)dx 


B 


li 


l n 

, O, 


[ d * w °w l i x - 1 1 

r ( dw°(x)Y 

J dx 2 2 J 

{ dx ) 


( 171 ) 


Note that the solution for u °(x) involves the three terms which are the integral of w°(x) , 
its first derivative, and the integral of the square of its first derivative. Since w°(x) is 
an even function of x, these three terms are odd functions of x. Therefore, since 
symmetry about x = 0 has been assumed, and constants are even functions, 


C„ = 0 . (172) 

Since the solutions for w°(x) comprise three different functional forms which depend 
on the magnitude of the compressive axial load N, so will the solution for u °(x ) . The 
solution for v°(x) was presented in the previous section in eq. (142) and (143), and 
involves the first derivative of the solution for w*(x), given by eqs. (166), (167), and 
(168). 

Since the expression for u °(x ) , eq. (171), involves the integration of the square of 
the first derivative of w®(jc) , the solution is not easily obtainable. However, after much 
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algebra, it can be shown that, 


for 0 zN>N* , 


u°(x) = CjX + C 2 sinh(ax)cos(P*) + C 3 cosh(ax)sin(px) 
+ C 4 sinh(2ax)cos(2px) + C 5 cosh(2ax)sin(2px) 

+ C 6 sinh(2ax) + C 7 sin(2Px) 


where 



\/ 


‘12 


^ 11^22 ^12 


n+n: 


» ) 


l 12 


^ 11^22 ^12 


-nI+ 


FGa p 
2 


B n (Fa+GP) + ^F(P-a) 

A 

A n {a 2 + P 2 ) 


C 2 = 


fl n (Ga-Fp) + ^F(a + P) 

c ’ *„(«»♦« 

(g 2 -f 2 )« p 2 -fg(« 2 p + p 3 ) 

4 " 8(a 2 + P 2 ) 

(G 2 -F 2 )g 2 p-FG(tt 3 + «p 2 ) 
3 8(a 2 + P 2 ) 

C 6 = - (F 2 a 2 + G 2 p 2 )/ (8 a) 

c (F 2 a 2 + G 2 P 2 )(« 2 + P 2 ) 

7 8(a 2 P + P 3 ) 


and 


v^x) = — | F[asinh(ax)cos(Px) - pcosh(ax)sin(Px)] 
^66 ^ 

- G[acosh(ax)sin(Px) - Psinh(ax)cos(Px)] 



(173a) 


(173b) 
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For N = AT, 


«°(x) = CjX + C 2 x 3 + C 3 sin( Px) + C 4 sin(2 px) 

+ C 5 xcos(Px) + C 6 xcos(2Px) + C 7 x 2 sin(2Px) 


where 


C.= 


,2 V 


1 + 


‘12 


^ 11^22 ^12 


n+n: 


V A n / 


^12 n t^ F 2 p 2 +FGp-G 2 


^ 11^22 ^12 
^ ' 


2 o2 


G Z P 


12 


C 3 = 




F- — | 

L_L 2 


a ll 


c _ (2G 2 -2F 2 p 2 -FG P) 
4 ~ 16 p 


in rn + ^ l2 


mi 


C 6 = -FGp/8 


c 7 = -G 2 p/8 


and 


v°(x) = 


^ j F [-Psin(px)j + G [sin(Px) + pxcos(Px)] } 
A 66 1 > 


S 

+ X 

*66 


(174a) 


(174b) 
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For N < N\ 


u °(x) = C x x + C 2 siiijp,*) + CjSinfP^x) + C 4 shi(2P I x) + C^in^P^) 
+ C 6 sin[(pj - P 2 )*] + C 7 sin[(P, + P 2 )x] 


where 


Ci- 


1 +- 


l 12 


V iV+M n 


^ 11^22 ”^12 y 


V 

< A n 

-BnFVr 


\ A n ) 


A 12 .,r J*P? + G*P* 


12 Ne - 


C 

U 2 


/? P 1 


^ 11^22 ~^12 


c," 


4 

* p 2 


*11 


*11 


C 4 = 


C 6 = - 


F 2 Pi 

8 

FGPA 

4(Pi-P 2 ) 


G 2 P 2 


C ’ = 8 


- ^P.P2 

4(Pi-P 2 ) 


(175a) 


and 


(*) ■ ^{fl-p.sii^p,*)] * G [-PjSii^pjX)] | + j-J 
■^66 ^ 


(175b) 


The axial compressive load N * , which corresponds to the load at which the character 
of the roots to the characteristic equation, eq. (152), changes from the roots being 
complex conjugates to the roots being repeating pure imaginary roots, has been shown 
to correspond to the collapse load of the cylinder by Booton (ref. 2). Since 
application of a compressive axial load corresponding to the collapse load would cause 
catastrophic failure of the cylinder, the analyses and results discussed in the remainder 
of this work will focus only on compressive axial loads in the range 0 i N> N* . 
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E. Preloading Response Due to Thermal Effects 

Composite cylinders are usually fabricated on a male mandrel and are consolidated 
at an elevated temperature. After consolidation, the temperature is lowered to ambient 
temperature and the male mandrel is removed from the cured cylinder. If the fabrication 
and consolidation are assumed to be axisymmetric, the cured shape of the cylinder can 
be determined using the solutions from the previous sections. 

Since there are no loads applied to the cylinder after the mandrel is removed, N is 
set equal to zero, as is S, and the response is given by eq. (166). Also, the roots are 
given by eqs. (158) and (157). Note that the roots are only functions of the material 
properties and cylinder geometry, since N is set equal to zero. Therefore, the shape of 
the deformed cylinder is a function of the material properties and geometry only. The 
particular solution, eq. (150), which is a function of the material properties, cylinder 
geometry, and the temperature change, governs the radial deformation of the cylinder 
away from the ends. 

Since the boundaries of the cylinder are unrestrained after the mandrel is removed, 
the appropriate boundary conditions for this case are the lubricated boundary conditions, 
eq. (146), with N set equal to zero. 

In summary, from the material properties and geometry of the cylinder, the root parts 
a and P of eq. (166) are known. By enforcing the lubricated end boundary condition, 
the constants F and G of eq. (166) can be solved using eqs. (146). From the material 
properties, cylinder geometry, and temperature change, the particular solution, eq. (150), 
with N set equal to zero, is known. Therefore, the deformed shape w°(x) can be 
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calculated. 


F. Numerical Results for the Case of Thermally-Induced Preloading 

Relations between - — , — ~ ff — and j for three 16 layer cylinders are 

H H H L, 

presented in Fig. 11, Fig. 12, and Fig. 13, respectively. These cylinders have stacking 
sequences of [+45/-45/0 ? k<;. [+45/-45/0J4t> and [02/ _ 45/+45] 4T , a length to radius 
ratio, L/R, of 3 and radius to thickness ratio, R/H, of 125. The value of these 
parameters are representative of thin, moderately long cylinders. The results are felt to 
be valid for any cylinder with L/R ^ 2 and R/H ^ 100. Specific dimensions used 
herein can be determined knowing a single layer of fiber reinforced material is 0.005 in. 
thick. The layer material properties used in the calculations are given in Table I. 


Table I. Layer Material Properties 


Ei 

(Msi) 

E* 

(Msi) 

Gij 

(Msi) 

*12 

layer 
thickness, 
h (in.) 

«i 

(in./in.)/°F 

<*2 

(in./in.)/°F 

20. 

1.3 

1.03 

.3 

.005 

-.167 x 10* 

15.6 x 10* 


The first stacking sequence represents an often-used slightly orthotropic lay-up, while 
the second and third stacking sequences represent two unsymmetric variants of the first 
one. If unsymmetric laminates are to gain favor, slight deviations from symmetry, as 
with the second lam inates, are most likely to initially be used. The laminate properties 
and thermally-induced equivalent stress resultants with AT=-280°F for these three 
cylinders are given in Table n. 
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Table II. Laminate Properties and Thermally-Induced Stress Resultants for 
[+45/-45/0J2S, [+45/-45/0J 4t , and [02/-45/+45] 4T Cylinders, 

AT=-280°F. 


Laminate Property 

[+45/-45/02] 2s 

[+ 45 /- 45 /O 2 LT 

[0 2 /-45/+45] 4T 

A„ (lb/in) 

1.068 x 10 6 

1.068 x 10 6 

1.068 x 10 6 

Aj 2 (lb/in) 

.1966 x 10 6 

.1966 x 10 6 

.1966 x 10 6 

A 22 (lb/in) 

.3156 x 10 6 

.3156 x 10 6 

.3156 x 10 6 

A* (lb/in) 

.2476 x 10 6 

.2476 x 10 6 

.2476 x 10 6 

B„ (lb-in/in) 

0 

2,707. 

- 2707. 

B ,2 (lb-in/in) 

0 

-826.0 

826.0 

B 16 (lb-in/in) 

0 

-470.3 

470.3 

D n (lb-in 2 /in) 

461.3 

569.6 

569.6 

N X T Gb/in) 

-160.2 

-160.2 

-160.2 

N/ Gb/in) 

-357.1 

-357.1 

-357.1 

M X T Gb-in/in) 

0 

.4921 

-.4921 


O [+45/-45/02I2S 

0.0100 - 0 [+45/-45/O2UT 

A (02/-45/+45UT 


0.0025 


0.2 0.3 

x/L 


Fig. 11. Dimensionless Axial Mid-surface Displacement of Cylinders with 
Lubricated Boundary Conditions, N=0, AT=-280°F. 
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Fig. 12. 


Fig. 13. 
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Dimensionless Tangential Mid- surface Displacement of Cylinders 
with Lubricated Boundary Conditions, N=0, AT=-280°R 
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Dimensionless Radial Mid-surface Displacement of Cylinders with 
Lubricated Boundary Conditions, N=0, AT =-280 F. 


Fig. 11 reveals that the axial displacement of the cylinder mid-surface for all three 
cylinders is essentially a linear function of x and that the magnitudes of the responses for 
these cylinders is the same. Since the effective axial coefficient of thermal expansion of 
the cylinders is negative, they respond to the temperature change of AT=-280°F by 
expanding axially. 

By examination of Fig. 12 and Fig. 13, it is seen that for the symmetric cylinder, the 
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tangential displacements of the mid-surface are zero over the entire cylinder length, and 
radial displacements of the mid-surface of the symmetric cylinder are constant. 

Comparing the results of the unsymmetric [+45/-45/0J 4T cylinder to the results of 
the [(y-45/ +451 , t cylinder, it is seen that the tangential and radial displacements vary 
along the half-length of the cylinders, particularly toward the cylinders’ ends. The ends 
of the unsymmetric cylinders curl radially and twist, producing a boundary layer effect. 
The direction of the radial curl is a function of the sign of the thermally induced moment 
Af/, resulting from the opposite stacking sequences of the cylinders. The [+45/-45/0J 4T 
cylinder curls outward, while the [0 2 /-45/+45] 4T cylinder curls inward. Fig. 12 
illustrates that both unsymmetric cylinders twist in the same direction due to the 
thermally-induced preloading, although the [0 2 /-45/+45] 4T cylinder exhibits a slightly 
larger tangential displacement at the end. It is important to realize that the study of 
unsymmetrically laminated cylinders would not be correct without including this 
thermally-induced deformation due to cooling from the consolidation temperature. 

G. Cylinder Response Due to Thermally-Induced Preloading Effects and a 

Compressive Axial Load 

In order to correctly model the response of the cylinder due to thermal preloading 
effects and a compressive axial end load, the thermally-induced preloading deformation 
at the end of the cylinder must be taken into account when enforcing the end conditions 
under axial load. Since the cylinder deforms due to cooling from consolidation 
temperature to ambient temperature, any fixture used to apply the axial end load to the 
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thermally deformed cylinder must take into account the thermally-induced tangential and 
radial deformation and slope at the end of the cylinder before loading. This impacts the 
specification of the simply supported and clamped boundary conditions, eq. (147) and 
(148), respectively. Stated in another way, simply supported or clamped end conditions 
would resist any tangential and radial end displacements relative to the thermally-induced 
preloading value, and clamped end conditions would also resist rotation of the ends 
relative to the thermally-induced preloading value. As observed in the previous section, 
under thermally-induced preloading effects, for unsymmetric laminates, the radius at the 
end of the cylinder will most likely be different from the radius of the cylinder away 
from the ends. Therefore, the axial load N would be applied eccentrically relative to 
the mid-length of the cylinder. This could have an influence on the response of the 
cylinder near the ends. 

Since the boundary conditions associated with the thermal-induced preloading effects 
are those of lubricated ends, it follows that the thermally-induced torsional load, S T , is 
zero. Therefore, the tangential displacement v° at the end of the cylinder after 
thermally-induced preloading effects is given by eq. (142) with S equal to zero, i.e., 



B^dw£\ 

*66 


( 176 ) 


Under ther mall y-induced preloading effects and axial load, the simply supported and 
clamped boundary conditions are: 
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1) Simply Supported end (at x = +L/2): 
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(177c) 


2) Clamped end (at x=+L/2): 
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and 


dw' 


dx 


dwj 




dx 


-I 


(178c) 


where the subscript T denotes the value of the response due to thermally-induced 
preloading effects. 

Obviously, since w*(x) and its first derivative are involved in the boundary 
conditions, the unknown constants F and G in the solutions for w°(x), eqs. (166), 
(167), and (168), can be solved for. Also, these two boundary conditions require that 
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a third unknown, the torsional load S be found. Therefore, three unknowns, F, G, and 

5, must be solved for using eqs. (177) and (178). 

It is worth noting that in the case of clamped boundary conditions, the specification 

. o dwj 

of v®, eq. (178a), and specification of — — - , eq. (178c), both involve —j— evaluated 

« B 

at x = +L/2. In fact, if eq. (178c) is multiplied by and the result is subtracted from 

A 66 

eq. (178a), the solution 

5 = 0 < 179 > 

is found. Therefore, eq. (178a) can now be eliminated from the system of equations 

used to solve for F and G, leaving eq. (178b) and (178c) as the system of equations to 

be solved. The physical interpretation of this result is that for these particular 

unsymmetrically laminated cylinders (which result in certain B matrix terms not being 

zero), axial compression with clamped boundary conditions does not induce a torsional 

load 5 , while simple support boundary conditions will induce a nonzero S . This is a 

dw Q 

result of the definition of the tangential displacement v®(x) involving the slope and 
the value of the thermally-induced preloading torsional load, S T , being equal to zero. 
If the boundary conditions during the cooling from consolidation temperature to ambient 
temperature were other than lubricated end conditions, i.e., if the ends were simply 
supported or clamped during the cooling procedure, a nonzero value for S T would result 
for unsymmetrically laminated cylinders, and imposition of clamped boundary conditions 
during axial compression would result in an induced torsional load 5. 
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H. Numerical Results for the Case of Thermally-Induced Preloading Effects and a 
Compressive Axial Load 

Dimensionless axial, tangential, and radial displacements for cylinders with 
[+45/-45/0J2S, [+45/-45/0J4T, and [(y-45/ +45] 4T stacking sequences, with simply 
supported and clamped end conditions and varying load N, are presented in Fig. 14 
through Fig. 31. The figures illustrate the variation of these displacements along the 
dimensionless half-length of the cylinders. The three families of figures illustrating the 
axial, tangential, and radial mid-surface displacement have common vertical scales for 
ease of comparison within these families. 

The two axial load levels investigated are given by fractions of the load N * , i.e., 
for N= 10% N* and N= 90% N * . The quantity N* is independent of boundary conditions 
and thermal preloading. Recall that the quantity N* dictated the form of the roots of the 
characteristic equation, eq. (152). The values of N* for the three cylinders are: 

[+45/-45/0J2S : -2271 lb/in, 

[+45/-45/0J4T : -2241 lb/in, 

[<V-45/+45] 4T : -2771 lb/in. 

L Simolv Supported Boundary Conditions with Thermallv-Induced Preloading Effects 
and a Compressive Axial Load 

The dimensionless mid-surface displacements as a function of distance along the half 
length of these cylinders with simply supported aids are presented in Fig. 14 through 
Fig. 16 for N=10% N * , and Fig. 17 through Fig. 19 for N= 90% N * . Comparison of 
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the results relative to the two load levels exhibits the change in the magnitude and shape 
of the responses. The difference in the shape of the responses is a result of the nonlinear 
nature of the governing equations. In particular, the length of the boundary layer 

increases with increasingly compressive axial load. 

Fig. 19 illustrates that, at comparable load levels, the maximum value of the radial 
response of the [+457-45/0^ cylinder is significantly larger than that of the 
[(y-45/+45] 4T cylinder. Note that the actual load N = 90% N* on the [+45/-45/0J4T 
cylinder has a smaller magnitude than the load on the [02/-45/+45] 4T cylinder, since the 
load N* is larger in magnitude for the [0 2 /-45/+45] 4T cylinder and each cylinder in these 
figures is loaded axially based on the quantity N/N * , with N" being different for each 
cylinder. 

The tangential displacement v°(x) is given by eq. (142). Note that the magnitude of 

i4 w and the magnitude of B l6 (presented in Table II) for the [+45/-45/0J 4T and 

[ 02 /- 45 /+ 45 ] 4T cylinders are identical, while the sign of B 16 for these two cylinders are 

opposite. Therefore, the difference in the magnitude of the tangential displacement, 

ctw Q 

presented in Fig. 15 and Fig. 18, is a result of the difference in the magnitude of 
for these cylinders, while the difference in the sign of v*(x) is a result of the difference 
in the sign of B l6 for these cylinders. 

As illustrated in Fig. 14 and Fig. 17, the magnitudes of the axial displacements of 
the [(V- 45 /+ 45 ] 4 t cylinder are larger than the magnitudes of the axial displacements of 
the [+45/-45/0J 4T and [+45/-45/0J2S cylinders. Since all three cylinders have identical 
inplane stiffnesses, as dictated by A u , and and the fact that these figures 
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represent displacement data based on dimensionless loads N/N' , this difference in 
magnitude is a consequence of the larger magnitude of N * for the [(y- 45 /+ 45 ] 4 T 
cylinder. 



Fig. 14. Dimensionless Axial Mid-surface Displacement of Cylinders with 
Simply Supported Ends, N=10%N*, AT=-280°F. 
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Dimensionless Tangential Mid-surface Displacement of Cylinders 
with Simply Supported Ends, N=10%N*, AT=-280°F. 
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Fig. 16. Dimensionless Radial Mid-surface Displacement of Cylinders with 
Simply Supported Ends, N=10%N\ AT=-280°F. 
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2, Clamped Boundary Conditions with Thermally-Induced Preloading Effects and a 

Compressive Axial Load 

Since simply supported boundary conditions are difficult, if not impractical, to obtain 
in reality, and since the displacements of unsymmetrically laminated cylinders are large 
and rapidly changing near the ends of the cylinder, it is legitimate to ask if the responses 
of cylinders with clamped ends would be comparatively diminished. 

The dimensionless mid-surface displacements as a function of distance along the half 
length of cylinders with clamped ends are presented in Fig. 20 through Fig. 22 for 
N=10% AT, and Fig. 23 through Fig. 25 for N=90% AT. 

As was the case with simply supported ends, it is observed that the axial displacement 
is largest for the [02/-45/+45] 4T cylinder. Again, this can be attributed to the fact that 
the value of the axial load N is significantly larger for this cylinder, relative to the other 
two, while the inplane stiffnesses of all three cylinders are identical. 

From Fig. 24 and Fig. 25, it is apparent that the maximum values of the tangential 
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and radial displacements of the [0 2 /-45/+45] 4T cylinder are larger than those of the 
[+45/-45/0J4T cylinder. This is converse to the results of the previous section, Fig. 18 
and Fig. 19, where simply supported boundary conditions were imposed. In comparing 
the simply supported and clamped boundary condition results for N = 90% N , it is also 
evident that the range of values for the tangential and radial response of the 
[+45/-45/0J4T and [(y-45/+45]4 T cylinders with simply supported boundary conditions 
are broader than the range of values of these responses for the clamped boundary 
conditions. 


Again, from Fig. 22 and Fig. 25, it can be observed that the length of the boundary 
layer increases with increasingly compressive axial load. 
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Fig. 21. 



Dimensionless Tangential Mid-surface Displacement of Cylinders 
with Clamped Ends, N=10%N*, AT=-280°F. 



Fig. 22. Dimensionless Radial Mid-surface Displacement of Cylinders with 
Clamped Ends, N=10%N*, AT=-280°F. 
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Fig. 23. Dimensionless Axial Mid-surface Displacement of Cylinders with 
Clamped Ends, N=90%N*, AT=-280°F. 
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Fig 24. Dimensionless Tangential Mid-surface Displacement of Cylinders 
with Clamped Ends, N=90%N‘, AT=-280°F. 



^ Effect of Neglecting to Inclu d e the Thermallv-Induccd PreR>adWg Effects 

In Fig. 26 through Fig. 31, the dimensionless axial, tangential, and radial mid-surface 
displacements are presented for simply supported cylinders, subjected to the same axial 
loads as before, but neglecting to include the thermally-induced preloading effects. 

Comparing the case of simply supported cylinders with thermal preloading and axial 
load to the case of axial load only for N=90%N\ i.e.. Fig. 17 through Fig. 19 and 
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Fig. 29 through Fig. 31, it is evident that neglecting to include the thermally-induced 
preloading effects has a measurable effect on the predicted radial displacement response 
of the cylinders. There is not much difference in the axial or tangential displacement 
response. The tangential and radial responses at the mid-length of the cylinder differ 
between these two cases since the particular solution depends on the thermally-induced 
stress resultants Nj, n£, and Mj. 

Fig. 32 and Fig. 33 represent the displacement responses for a [Og^gJr cylinder 
using the same material properties as above. These figures illustrate the exaggerated 
effect of neglecting to include the thermally-induced preloading effects in the solution for 
the axially loaded case. There is a large difference between the results. This particular 
case was studied in ref. 1 without including thermal effects. Results such as shown in 
Fig. 32 and Fig. 33 certainly provide motivation to recompute the results of ref. 1 with 
thermal effects. 



x/L 

Fig. 26. Dimensionless Axial Mid-surface Displacement of Cylinders with 
Simply Supported Ends, N=10%N\ AT=0°F. 
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Fig. 30. Dimensionless Tangential Mid-surface Displacement of Cylinders 
with Simply Supported Ends, N=90%N*, AT=0°F. 
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Simply Supported Ends, N=90%N\ AT=0°F. 
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Fig. 32. Dimensionless Radial Mid-surface Displacement of a [0 g /90Jr 
Cylinder with Simply Supported Ends, N=10%N\ AT=0°F and 
AT=-280°F. 
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Fig. 33. Dimensionless Radial Mid-surface Displacement of a [0 8 /90 8 ]t 
C ylinder with Simply Supported Ends, N=90%N\ AT=0°F and 
AT=-280°F. 

In the next chapter, the effects of unsymmetric stacking sequence and thermally- 
induced preloading on the intralaminar stress components will be investigated. The 
solutions to the axisymmetric problem presented thus far will be employed to compute 
the intralaminar stress components in the individual layers of cylinders. 
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IV. CALCULATION OF INTRALAMINAR STRESSES 


In this chapter the equations used to calculate the principal material system stresses 
o u , Ojj , and t l2 within the wall of the [+45/-45/0J2S , [+45/-45/0J4T , and 
[0 2 / -45/ +45 ] 4t cylinders will be presented, and numerical results with thermally-induced 
preloading effects and compressive axial load will be illustrated. The principal material 
system stress components o„ , , and t 12 will be referred to as the principal material 

system intralaminar stresses. These stress components can be calculated from the 
principal material system strains, e n , e^, and Yj 2 » through constitutive relations. The 
other three stress components, t^, t 6 ,, and o r , are referred to as interlaminar stresses 
and will be discussed in a later chapter. 

A. Equations Describing Intralaminar Stresses 

Under the axisymmetric assumption of the previous chapter, the intralaminar strains 
t x , e e , and in the cylinder coordinate system x-0-r are given by eq. (12) and eq. 


(128), namely, 



and 


P'e'° 
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\ n o\ _o w‘ o 
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(128) 


In the previous chapter closed-form solutions for k°(x), v *(*) > w (*) were 
presented. Therefore, the derivatives of these functions with respect to x are analytically 
obtainable, and the expressions in eq. (128) are known for a given temperature change, 


axial load, and boundary conditions. 

The principal material strains can be calculated from the strains t x , e e , and Y*e by 
coordinate transformation from the x-0-z coordinate system to the 1-2-3 principal material 
coordinate system for each layer. These transformations are given by 
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(180) 


where 8 is the angle measured from the x axis to the fiber direction (the 1-axis) of a 


given layer. 

Once the total strains in the 1-2-3 coordinate system are known, the mechanical 
strains are obtained by subtracting the free thermal strains, 

e n ■ «n Ar » (181) 
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from the total strains given by eq. (180). By definition, an orthotropic layer does not 
have a free thermal shearing strain. Therefore, the mechanical strains are given by 

®n = ®n “ a n AT , 

4 = ® 22-«22 (182) 

M 

Y 12 = Y 12 • 

From these mechanical strains, the principal material system intralaminar stresses 
o a , Ojj , and Tjj can be calculated using the reduced stiffness matrix for each layer, 
i.e., 
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The reduced stiffnesses are 


(183) 


<?n = 

<?i2 = 


1 - v u v 2l 
V 12^2 

1 ~ V 12 V 21 ’ 

<?*- . 

* V 12 V 21 

= ®11 » 


(184) 


where E t , E 2 , v u , and G n are known for a given orthotropic layer and v 21 is 
calculated by 
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The values for the layer material properties used in these analyses were given in the 
previous chapter in Table I. 

B. Numerical Results for Intralaminar Stresses: Case of Thermally-Induced 

Preloading Effects and a Compressive Axial Load with Clamped Boundary 

Conditions 

Since clamped boundary conditions more closely represent actual applications of 
cylinders, the remainder of the analyses will focus on clamped boundary conditions. 
Again, the three cylinders of interest have a length to radius ratio of 3, a radius to 
thickness ratio of 125, and have [+45/-45/0J2S » [+45/-45/0J4T * and [02/-45/+45] 4T 
stacking sequences. The thermal preloading effects are due to a temperature change, 
AT, equal to -280°F. 

The numerical results to be presented in this chapter are principal material system 
intralaminar stresses in various groups of layers within the cylinder wall. The stresses 
are normalized by the quantity (N/H), where N is the compressive axial load and H 
is the wall thickness, 0.080 in. The quantity (N/H) represents the average axial stress 
in the cylinder wall. The stress components are reported at the mid-thickness location 
of each layer. The two axial load levels investigated are given by fractions of the load 
N\ i.e., N= 10% AT and N = 90%AT. 

For each cylinder and load case, six figures are presented. There are two figures for 


each of the mid-layer principal material system intralaminar stress components, o n , a n , 
and t 12 ; one figure for the 0° layers and one figure for the ±45° layers. Due to the 
difference in magnitude of the respective stress components in these two groups of 
layers, two figures can be used. For each cylinder and load case the vertical scale of 
each of the six respective figures is the same for comparative purposes. Since the 
magnitudes of the stress components are largest near the ends of the cylinders, the 
horizontal scale includes only the portion of the half-length of the cylinder from x/L=0.3 
to x/L=0.5. Where possible, the layer number in the stacking sequence is printed next 
to the respective layer’s stress component relation in the figures. Where the relations are 
too closely spaced in the figures, the direction of increasing layer numbering in the 
stacking sequence is shown with arrows. Recall, layer no. 1 is at the inner most radial 
position. 

The first three sets of figures, Fig. 34 through Fig. 36, represent principal material 
system intralaminar stress results for the three cylinders subjected to a compressive axial 
load of N ■ 10% N\ including thermally-induced preloading effects, i.e., with a 
temperature change, AT, equal to -280°F. The next three sets of figures, Fig. 37 
through Fig. 39, represent principal intralaminar stress results for the three cylinders 
subjected to a compressive axial load of N - 90% W*, including thermally-induced 
preloading effects. The last set of figures, Fig. 40, represent principal intralaminar stress 
results for the [+45/-45/0J4T cylinder subjected to a compressive axial load of 
N = 90%N* , neglecting thermally-induced preloading effects, i.e., with a temperature 
change, AT, equal to 0°. 
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C. Discussion of Intralaminar Stress Results 

For the low load level for the [+45/-45/0J2S laminate, Fig. 34, the fiber-direction 
stresses, o u , are higher in the ±45° layers than the 0° layers. Since the 0° layers are 
aligned with the applied load and thus will directly bear the load, it is surprising that the 
fiber-direction stresses are not higher in the 0° layers. This anomaly is no doubt due to 
the thermally-induced effects dominating at this low load level. This will be seen 
shortly. The intralaminar stresses perpendicular to the fibers, o^, are about the same 
for the 0° layers as for the ±45° layers. The shear stresses are small, being zero in the 
0° layers and equal and opposite in the ±45° layers. 

For the low load level, the stresses in the [+45/-45/0j4 T laminate, Fig. 35, are 
similar in magnitude and spatial distribution to the stresses in the symmetric laminate. 
The primary difference between the unsymmetric laminate and the symmetric laminate 
is a minor level of shear stress in the 0° layers of the unsymmetric laminate. 

The stress levels at the low load level in the [(y-45/+45] 4T laminate, Fig. 36, are 
similar to the stresses in the [+ 45 /- 45 /O 2 I 2 S and [+45/-45/0j4 T laminates. The fiber- 
direction stresses in the 45° layers, at the end of the cylinder, are somewhat lower for 
the [(y-45/+45]4T laminate than for the other two laminates. The compressive stresses 
perpendicular to the fibers, , are not as high as for the [O 2 /- 45 / H~ 45) 4T laminate as for 
the other two. The same is true for the shear stresses in the ±45° layers, though the 
difference in stress levels for one laminate to the other two is of little consequence for 
this low load level. 

For the high load level and the symmetrically laminated cylinder. Fig. 37, the fiber- 
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direction stress levels in the 0° layers, the load bearing layers, are higher than in the 
±45° layers. Apparently, at the low load level, Fig. 34a and Fig. 34b, thermal effects 
do dominate but at the higher load level, they are in the background. This will be 
discussed more shortly. 

The stress levels in the two unsymmetric laminates at the high load are quite similar 
to the stress levels in the symmetric laminate at that load level. The intralaminar stresses 
perpendicular to the fiber-direction, a n , are compressive, thus virtually eliminating the 
potential for micro-cracking due to this stress component. The shear stresses in the 
±45° layers for the three laminates could lead to matrix cracking, but they are low, even 
for this high load level. 

To provide an indication of the magnitude of the thermally-induced intralaminar 
stresses, the stresses in the [+45/-45/0J 4T cylinder with the high load level but not 
including thermal effects are shown in Fig. 40. These should be compared directly with 
Fig. 38, the same laminate and same load level, but with thermal effects included. All 
stresses except the fiber-direction stress in the 0° layers shows some influence of the 
thermal effects. By comparing Fig. 38c and Fig. 38d with Fig. 40c and Fig. 40d, it 
appears as though the intralaminar stresses perpendicular to the fibers, o^, due to 
thermal effects, are compressive. Also, when comparing Fig. 38f with Fig. 40f, it can 
be concluded that thermal effects relieve somewhat the intralaminar stress in the ±45° 
layers. 

It should be noted that the solutions for the intralaminar stresses obtained to this point 
can be used directly in a failure criterion associated with classical lamination theory 
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(CLT). However, since the failure of a cylinder on a material failure basis also consists 
of determining whether delamination occurs for a given load case, a method for 

determining the interlaminar stresses is needed. 

In the following chapters, the analyses performed thus far will be used to derive 
solutions for the interlaminar stress components within the cylinder wall. 
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Fig. 34a. a n along x/L for 0° Layers of 
a [+45/-45/0J2S Cylinder, N=10%N*, 
AT=-280°F. 
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Fig. 34b. (r n along x/L for 45° Layers 
of a [+45/-45/0J2S Cylinder, N= 10%N*, 
AT=-280°F. 
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Fig, 34c. 022 along x/L for 0° Layers of 
a [+45/-45/02J2S Cylinder, N=10%N‘, 
AT=-280°F. 


Fig. 34d. a 2 2 along x/L for 45° Layers 
of a [+45/-45/0J2S Cylinder, N=10%N*, 
AT =-280°F. 
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Fig. 34e. r, 2 along x/L for 0° Layers of 
a [+45/-45/0J2S Cylinder, N=10%N\ 
AT=-280°F. 
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Fig. 34f. r n along x/L for 45° Layers of 
a [+45/-45/02]2s Cylinder, N=10%N*, 
AT=-280°F. 
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Fig. 35e. r 12 along x/L for 0° Layers of 
a [+45/-45/0J4T Cylinder, N=10%N\ 
AT=-280°F. 



Fig. 35b. <r„ along x/L for 45° Layers 
of a [+45/-45/0J 4T Cylinder, N=10%N*, 
AT=-280°F. 



Fig. 35d. on along x/L for 45 ° Lay 6 ™ 
of a [+45/-45/0J4T Cylinder, N=10%N*, 
AT=-280°F. 
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Fig. 35f. r, 2 along x/L for 45° Layers of 
a [+45/-45/0J4T Cylinder, N=10%N*, 
AT=-280°F. 
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Fig. 36a. <r„ along x/L for 0° Layers of 
a [02/-45/+45]4 T Cylinder, N=10%N’, 
AT=-280°F. 
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Fig. 36b. a n along x/L for 45° Layers 
of a [(y-45/+45]4 T Cylinder, N = 10%N\ 
AT=-280°F. 
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Fig. 3oCe 022 along x/L for 0° Layers of 
a [(y-45/+45] 4T Cylinder, N=10%N*, 
AT=-280°F. 


Fig. 36d. a n along x/L for 45° Layers 
of a [(V-45/+45]4 T Cylinder, N= 10%N*, 
AT =-280°F. 
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Fig. 36e. Tj 2 along x/L for 0° Layers of 
a [(V-45/+45LT Cylinder, N=10%N*, 
AT=-280°F. 
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Fig. 36f. r 12 along x/L for 45° Layers of 
a [(V-45/+45LT Cylinder, N=10%N*, 
AT=-280°F. 
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Fig. 37a. <r„ along x/L for 0° Layers of 
a [+45/-45/0J2S Cylinder, N=90%N*, 
AT=-280°F. 



Fig. 37c. a-n along x/L for 0° Layers of 
a [+45/-45/0J2S Cylinder, N=90%N*, 
AT=-280°F. 
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Fig. 37e. Tj 2 along x/L for 0° Layers of 
a [+45/-45/0J2S Cylinder, N=90%N*, 
AT=-280°F. 
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Fig. 37d. o-a along x/L for 45° Layers 
of a [+45/-45/0J2S Cylinder, N=90%N\ 
AT=-280°F. 
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Fig- 37f. r 12 along x/L for 45° Layers of 
a [+45/-45/0J2S CyUnder, N=90%N*, 
AT=-280°F. 
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Fig. 38e. Tj 2 along x/L for 0° Layers of 
a [+45/-45/0J4T Cylinder, N=90%N*, 
AT =-280°F. 


Fig. 38f. Tjj along x/L for 45° Layers of 
a [+45/-45/021,t Cylinder, N=90%N\ 
AT=-280°F. 
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Fig. 39a. o u along x/L for 0° Layers of 
a [0 2 /-45/+45] 4T Cylinder, N=90%N*, 
AT=-280°F. 


Fig. 39b. cr,, along x/L for 45° Layers 
of a [(y-45/+45] 4T Cylinder, N=90%N\ 
AT=-280°F. 
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Fig. 39c. 022 along x/L for 0° Layers of 

a [(V-45/+45] 4T Cylinder, N=90%N*, 
AT=-280°F. 
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Fig. 39d. ajj along x/L for 45° Layers 
ofa[(y-45/+45] 4T Cylinder, N=90%N\ 
AT=-280°F. 
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Fig. 39e. r 12 along x/L for 0° Layers of 
a [02/-45/+45] 4T Cylinder, N=90%N\ 
AT=-280°F. 


Fig. 39f. r, 2 along x/L for 45° Layers of 
a [0 2 /-45/+45] 4T Cylinder, N=90%N\ 
AT=-280°F. 
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Fig. 40a. <x n along x/L for 0° Layers of 
a [+45/-45/0J 4T Cylinder, N=90%N\ 
AT=0°F. 


Fig. 40b. a u along x/L for 45° Layers 
of a [+45/-45/0J4T Cylinder, N=90%N*, 
AT=0°F. 
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Fig. 40c. <T 22 along x/L for 0° Layers of 
a [+45/-45/0J4T Cylinder, N=90%N*, 
AT =0°F. 
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Fig. 40d. 0 j 2 along x/L for 45° Layers 
of a [+45/-45/0J4T Cylinder, N=90%N*, 
AT =0°F. 
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Fig. 40e. tj 2 along x/L for 0° Layers of 
a [+45/-45/0 J 4t Cylinder, N=90%N\ 
AT=0°F. 
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Fig. 40f. r u along x/L for 45° Layers of 
a [+45/-45/02 ] 4t Cylinder, N=90%N\ 
AT =0°F. 












V. DERIVATION OF THE THREE-DIMENSIONAL 
EQUILIBRIUM EQUATIONS IN CYLINDRICAL COORDINATES 

In the previous chapters, the assumptions of thin shell theory and axisymmetry have 
been employed in order to obtain closed-form solutions for the response of composite 
cylinders. This method has provided relations for the displacements and intralaminar 
strains and stresses as functions of the axial and radial coordinates. 

As a practical matter, it is important not to ignore the existence of the stresses at the 
layer-to-layer interfaces of the cylinder. These stress components are the shear stress 
components and x^, and the normal stress o r . These three stresses will be referred 
to as interlaminar stresses. If acceptable levels for the interlaminar stresses are 
exceeded, delamination of the layers can result, leading to a structural failure. The 
theory presented thus far incorporates classical lamination theory (CLT) which assumes 
a state of plane stress, i.e., that the interlaminar stresses can be ignored. Therefore, a 
different set of governing equations is required in order to calculate the interlaminar 
stresses. 

The aim of the remainder of this work is to calculate the interlaminar stresses by 
using the solutions for the intralaminar stresses derived in the previous chapters. The 
procedure used to this end consists of deriving the three-dimensional equilibrium 
equations and boundary conditions for a linear elastic body undergoing large 
deformations, simplifying these equations in accord with the pertinent assumptions of the 
previous chapters, and solving these equations using the solutions previously obtained for 
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the intralaminar stresses. 


In the following section, the three-dimensional equilibrium equations and boundary 
conditions will be derived for the cylindrical coordinate system used thus far in the 
analyses. In a subsequent section, the displacement gradients and rotations will be 
derived for the cylindrical coordinate system. The equilibrium equations and boundary 
conditions will then be recast in a form which explicitly illustrates the terms which 
depend on the displacement gradients and rotations, then, the equations can be simplified 
in a rational manner based on the magnitude of the deformations and rotations at a point. 
In order to bring the three dimensional equilibrium equations, boundary conditions, and 
the displacement gradients and rotations into accord with the solutions for the 
intralaminar stresses, these relations will be simplified under the assumptions of 
axisymmetry. In addition, the displacement gradients and rotations will be simplified 
further under the assumption of Kirchhoff, i.e., that the displacements vary linearly 
through the thickness of each layer. The simplified relations for the displacement 
gradients and rotations for the three cylinders of the previous analyses will be used to 
quantify the relative magnitudes of the coefficients of the stress components appearing 
in the three-dimensional equilibrium equations. These results will be used to further 
simplify the equilibrium equations. Finally, in the next chapter, these simplified 
equilibrium equations will be solved for the interlaminar stresses and results for the three 
cylinders will be presented. The method of solution will incorporate the solutions of the 
previous chapter for the intralaminar stresses, in the manner described by Pagano 
(ref. 3). 
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A. Transformation of the Three-Dimensional Equations of a Linear Elastic Body 
in Rectangular Coordinates to Cylindrical Coordinates 

The method of derivation of the three-dimensional equilibrium equations in cylindrical 
coordinates will begin with the derivation of the tensor form of the equilibrium equations 
in rectangular coordinates. These equations will then be transposed into the cylindrical 

coordinate system by methods of tensor analysis. 

The total potential energy (ref. 4) for a linear elastic body in rectangular 

coordinates is 

ntu ) = fJjWi - f ‘“) dv - f/‘“‘ • <18S| 


If Green’s Strain tensor, 


*1) = h u ij + u M +u p.‘ u pj) • 


(187) 


is incorporated into eq. (186), and the first variation is taken with respect to the 
displacement tensor u., the Euler equations and variationally consistent natural boundary 

conditions are 

inV > (188a) 


and 


h 4 ',V ’ h 00 S * (188b) 

OR Uj is specifiedoo S 2 , 


where T y are the Kirchhoff stress tensor components, u } are the displacement vector 
components, fj are the components of the body force vector, ij are the tractions specified 


117 


on the surface S 2 of the total boundary S, and V is the volume containing the body. As 
essential boundary conditions, the displacements are assumed to be zero over the surface 
Sj = S - S 2 , i.e., 

Uj = 0 on S x . (188c) 

The tensor tj is defined as 

tj = n. , (189) 

which acts over the surface S 2 , where n i are the components of the surface normal 
vector. Therefore, from the definition of the principal of minimum total potential 
energy, equations (188) represent the equilibrium equations and variationally consistent 
boundary conditions governing the three-dimensional response of a linear elastic body 
undergoing large deformations in a rectangular coordinate system. It should be noted 
that these are tensor equations, therefore, they can be transposed into any specific 
curvilinear coordinates in Euclidean space through the rules of tensor analysis. Note 
also, that the covariant and contravariant components of the tensor quantities in these 
equations are identical, since the coordinate system is rectangular. In the following, 
these equations will be transposed into curvilinear coordinates, and the distinction 
between covariant and contravariant tensor components must be observed. 

Eringen (ref. 5) notes that a tensor equation in rectangular coordinates can be 
resolved into curvilinear coordinates by enforcing the following two rules: "(a) The 
partial differentiation symbol (,) must be replaced with covariant differentiation (|), and 
(b) The repeated indices must be on diagonal positions." Applying these rules to eqs. 
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(188a) and (188b) results in a different form of the tensor equations, i.e., 

(o» 1- o"u'|,)l, ♦/' - 0 i”V, < 19#a) 


and 


t j + t p u J | = t 




on Sj , 


(190b) 


where covariant differentiation of contravariant tensors is now involved. The tensor 

calculus involved in carrying out these differentiations depends on the base vectors g k , 

( k\ 

the metric tensor g„, and the Christoffel symbols of the second kind < ) related to the 
transformation from rectangular to cylindrical coordinates. For the coordinate system 
used so far in the analysis, i.e., the coordinate system presented in Fig. 1, the base 


vectors are 


i 1 “ - 

g 2 = (rcos0)i 2 - (rsin0)i 3 
g 3 = (sin0)i 2 + (cos0)t 3 


(191) 


where i lf i 2 , and i 3 are the unit base vectors of the global rectangular coordinate 
system. Therefore, the metric tensor is 


llfvll = 


1 0 0 
0 r 2 0 , 
0 0 1 


( 192 ) 


and the nonzero Christoffel symbols of the second kind are 


Carrying out the differentiation in eq. (190a), according to the rules of covariant 
differentiation results in a tensor equation involving only partial derivatives, namely, 



The underscore of the superscripts of g u signify that the summation is suspended. In 
the following, the components of the body force, f J , are assumed to be zero. 

By expanding the tensor equation, eq. (194), for j = l, substituting the values of the 
metric tensor and Christoffel symbols derived above, collecting terms, and simplifying, 
the first equilibrium equation is 

r [ tU ( 1+ “u) + t ‘ 2 «u + 

+ r [ x21 ( 1+ “u) + t22 «U + t23 “lJ.2 (195a) 

+ [ r ( t31 ( 1 + «u) + t32 “u + t 33 “ J ].3 = 0 • 
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Likewise, for j=2, 


+ 


+ 


t 1 * 1 * 2,1 + t 12 (r 2 + ii^ + ru 3 ) + x l3 i* w - -^t 13 « 2 

T 21 1*2,1 + T 22 ^ 2 + tt 2>2 + ru 3 ) + T 23 «23 - “**“2 

X 31 « 21 + T 32 (r 2 + U 2> 2 + r ** 3 ) + * 33 «w ■ 7 t33 “2 

+ rx 21 1< 3 1 + x^ru^-Uj) + t 23 (r + r« 3>3 ) = 0 , 


(195b) 


and for j=3, 


T U l* 


3,1 


+ r 


x 21 u 


3.1 




- -ix 21 «. 


2,1 


+ T U ^*3,2 " ““ 2 j + tl3 ( 1 +W 3^) ^ 

^ 22 (« 3 .- 7 “ 2 ) + ‘ c 23 ( 1 + “ 33)] 2 
+ X 32 ^-^) + X 33 ( 1 + U 33))] 3 

" ^{ r 1 + + ru *) " " 7 “2) 


(195c) 


0 . 


It is important to realize that these equations involve tensor components x y and **, and 
not the physical components of the stresses and displacements. The physical components 
of stress will be denoted by o iJ as given by 





or 



(196) 


where 


1 

fyi = JJ 

8 “ 


(197) 


for orthogonal coordinate systems. Since the cylindrical coordinate system in this 
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analysis is orthogonal, then by eqs. (192) and (197), 


gil = 1 = gll , and 



(198) 


The physical components of the displacement tensor u < are u t as given by 




( 199 ) 


Using the above relations to obtain the physical components of the stress and 
displacement tensors, and 

(a) substituting x for the superscript 1, $ for 2, and r for 3 in eqs. (195), 

(b) substituting u for Uj , v for u 2 , and w for u 3 , and 


(c) replacing the partial derivative notation according to 


* ) - 


dx 




a more standard form of the equilibrium equations can be obtained, namely, 



(200a) 
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(200b) 


and 


a 

r — 

dv x j. 
— o + 

(. 1 dv w) 

1 + + — 

T *e + ^ X xr 

dx 

ax 

v r ae r j 

dr 

d 

+ — 

iv . i 

(. 1 dv w) 

0 8 ® + — T® r 

ae 

cbc 1 

^ r ae r j 

dr 


a 

+ — r 

^x” . ( 

! * I*: . j^x*' * 

iV' 

dr 

[a> ( 

rde r) 

dr 


^V 8 + -v'lo 98 * fi * — 

dx r\dQ J i dr 


dw\_ 9r _ 


jx 9 ' = 0 , 


dx 

_a_ 

ae 


d 

+ — r 

dr 


— o XI + -vW® + (l + ^V r 

dx r\dO J l 9rj . 

— x x9 + —[— -vlo 68 + (l + -^1 

dx r^ae I l dr) 

+ U^_ v ) x *r J l+ dw 

dx r\a0 j V dr 

. —X* 9 - (l + — — + — ) O #0 - — T er = 0 

dx { rdd r) dr 


eri 


^jo" 


(200c) 


In the last step, symmetry of the stress tensor was used in order to emphasize that six 
stress components, not nine, are involved. The fact that component designation of the 
stresses are superscripts is an artifact of the contravariant tensor of eqs. (195). However, 
since the stress components in eqs. (200) are physical components, they are not tensor 
components, and therefore, the terms covariant and contravariant have no meaning in this 
context. Hence, the difference between a superscript and subscript notation for the stress 
components is superficial. 

In a similar manner to that used for the equilibrium equations, the natural boundary 
condition equation, eq. (190b), can be transposed into the cylindrical coordinate system. 
The result, in the physical components of stress and displacement, is 
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and 


and 


du\ 

*} 

du\ 

*r 

du\ 


xO 


l^Ve + 

— T J 

r 30 

dr 

e e + 

du i 

r 30 

T 

dr 

+ 


r 30 

dr 


e r 


n. 


n r = t x on S 2 


OR u is specified on S, , 



1 + I* 

w' 
+ 

|t x0 ♦ 

3x ( 

r 30 


dr J 

♦ [ 

I + i* 

w' 
+ 

| c ee + dv x *r 

[dx \ 

r 30 


I dr 

\*L x xr + ( 

i + 

+ 

T® r + — rt rr li 


r 30 

H 

v “ V j 

3r 


n. 


- ;e 


on S, 


OR v is specified on S, , 


(201a) 


(201b) 


■£*" * t(w * ?)’"]*- 


ax 


=?r 


on S~ 


(201c) 


OR w is specified on S- , 


where the bar over the surface tractions t x , t e . and t r , and surface normals n , n.. 
and n r , in eqs. (201) signify that these are the physical components of these quantities. 

In summary, the equilibrium equations for the three-dimensional stress state of a body 
in cylindrical coordinates undergoing large deformation are three coupled, nonlinear 
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partial differential equations. The boundary condition relations are also coupled, 
nonlinear equations. 


B. Derivation of the Displacement Gradients and Rotation Components of the Finite 

Strain Tensor in Cylindrical Coordinates 
The equilibrium equations and boundary conditions just presented will later be 
simplified by rewriting them such that the coefficients of the stress components are cast 
in a form which reflects the contributions of the displacement gradients and rotations, in 
accord with the method presented by Novozhilov (ref. 6). This method enables 
one to simplify the equilibrium equations in a rational manner by eliminating certain 
terms based on the relative magnitudes of the products of the displacement gradients 
and/or rotations and stress components. 

First, the displacement gradients and rotations will be derived for the cylindrical 
coordinate system used thus far in the analyses. Enngen (Ref. 5) gives the following 
definitions for the displacement gradients, or infinitesimal strains, and the rotations, 

respectively, as 


where 



'u 


= '“'I*)’ 


(202a) 



(202b) 


If eq. (202b) is substituted into eq. (202a), the symmetry of the Christoffel symbols of 
the second kind is recognized, and the result is expanded, the following expressions for 
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the displacement gradient and rotations result 



r u 0 ( u k,i u i.k ) 


(203) 


e ik = e ki for all k4 ; r lk = -r kl for k*l , 

r kl = 0 for k= 0 . 

These definitions can be expressed in terms of their physical components by using the 
relation 


e ij ~ 


ig 


i U . 

g e “ 





(204) 


The displacements appearing in eq. (203) may also be expressed in terms of their 
physical components, as given by eq. (199). Therefore, the physical components of the 
displacement gradients and rotations in terms of the physical components of the 
displacements, can be shown to be 
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These relations can be expressed in the terminology of the cylindrical coordinate system 
by substituting, as before, the definitions u=u l , v=u 2 , and w = u 3 , and writing the 

ac ) 

= 1 i, . 

00 

These substitutions result in the following form of the displacement gradient and rotation 


partial differentials using the definitions - ( ),i » qq ( ),2 * ^ ^ 


definitions: 


= 


* A V** V/ r 

+ — 


A l/w v rr 

2 e = — + — 


du 
dx ’ 


CD 

s J. 

r 

0u _ 
00 ~ 

dv 

dx 

“ ~^ r ex 

1 du 
r 30 

0V 

+ dx ’ 

2 r xr 

du 

dr 

dw 

a c 

= 

~ 2 r rx » 

0u + 

0W 

0 r 

V 

0V 
+ 


0W _ _ 

dr 

’ 

xr e r 

r 

0r 

r 

00 ” 

1 

w 

+ — 






r 00 

» 

r 







1 0W _ V 






dr 

7 00 " r ’ 






dw 







dr 








e m 


- Vr M. vrr ' 

2 = — ♦ — 

0r 




( 206 ) 


The displacement gradients and rotations can be combined in order to make 
substitution into the equilibrium equations more obvious. Useful combinations of the 
displacement gradients and rotations are: 
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Substituting these combinations into the equilibrium equations, eqs.(200), derived in 
the previous section results in 

r + e **)° SX + ( e * + r *) T ‘ # + ( e ~ + r ") xX '] 

+ ^R 1 + + ^ + r xe)° ee + ( e *r + r -) T °1 (208a) 

+ ^ r R 1 + e xxK + i e x6 + + {*xr +r xr) arr ] = 0 * 
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* ^[('«e ” r »e) T “* * (* +e e») < ’” * (*«r + r «r) T * r ] (208b) 

- | '■[(«„ -r je )T" ♦ (1 - (*. r ‘r„)o"] 

* + («er->-e,) oM * (‘ + «")**' = 0 • 

a, ’ d * (1 *«„K'] 

- ^[(«„ - ♦ (<„ - r.,)<- M * (1 * «„)*•'] (208c) 

* * (‘«,-'-.r) t " * (‘ +'r')°"\ 

- - (> + '«)° M - («* ♦'*)**' - 0 • 

In a similar manner, the natural boundary condition relations, eqs. (201), can be 
rewritten by substituting the combinations of displacement gradients and rotations of eq. 
(207), i.e., 


[(> *‘„K + (*..* r -»K e + ( e » +r ..) t "F* 

+ [(> * '«) T ' 6 * ( e xe* r rt) 0 " " ( < »r* r xr) t, ']"e 
♦ [(1 ♦ e„)t- r ■* (« rt ♦ r rt )T e ' * («„ * r„)o"]I, - t * on S, 

OR u is specified on S 2 


(209a) 


and 



[ e xB ~ r xf>)° ZX 

♦( 

1 + e ee) T * 6 + («er + r er)^1»* 


+ 


♦( 

l + «ee)° 6e + ( e er + r er) xfl, ]«e 

A _ _‘1 4BPi ^ A jh 

(209b) 

♦1 


♦{ 

1 + «ee) T + ( e er + r e,)° ] n r = * 011 S 2 



OR v is specified on S 2 , 
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and 


+ («0, - r 0r) T ' 6 + (1 + *x 

+ [(*„- r ,rK e + ( V ^ r ) 0 ® 6 + (* + «„)**'] «6 

+ [{ e xr ~ r xr ) T ” + («0r "'er)* 8 ' + l 1 + «rr) a "] "r = *' 011 S 2 

OR w is specified on S 2 


(209c) 


To summarize the results of this section, relations between the displacement gradients 
and rotations and the partial derivatives appearing in the equilibrium equations for a 
cylindrical body undergoing large deformations, i.e., eqs. (200), and the natural 
boundary conditions, eqs. (201), have been obtained. These results were substituted into 
the equilibrium equations and natural boundary conditions derived in the previous 
section. In a subsequent section, the equilibrium equations and natural boundary 
conditions will be simplified through a rational method of comparison of the relative 
magnitudes of the displacement gradients and rotations in the coefficients of the stresses 
in these equations. 

C. Specialization of the Displacement Gradients and Rotations Under the 
Assumption of Axisymmetry 

Since the results of the previous chapters will be used to calculate the portions of the 
three-dimensional equilibrium equations just presented, the relations for the displacement 
gradients and rotations will be specialized under the assumption of axisymmetric response 
and loading, so that they will be in accord with the previous analyses. 
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The assumption of axisymmetric response and loading is defined by 


30 -o (210) 

00 


( ) being any response quantity. By enforcing this definition, the displacement gradients 
and rotations of the previous section, eqs. (206), become 




A vr 

2 «,e - “ 


du 
dx ’ 

dv 

dx ’ 

2e 

zr dr dx 
w 

> 

r 

0V _ V 

dr~ r ’ 
dw 

dr 




2 e Qr ~ 


e_. = 


2 ^*6 " 


2r„ = 


2 r 6r 


= -2r 


8bc 


e* » 


= -2r 

ar ax 

v dv « _ 

— + — = - 2 r. 


8r 


rO 


(211) 


D. Simplification of the Axisymmetric Displacement Gradients and Rotations Under 

the Assumptions of Kirchhoff and Donnell 
In order to incorporate the analyses of the previous chapters, the axisymmetric 
displacement gradients and rotations must be simplified by utilizing the relations 
describing the displacements from the previous chapters, i.e., 

u(x,e^) =«°(x,0) +zp"(x,6) 
v(x,e,r) = v °(x,6) + z Pe(x,6) 
w(*,6/) = w‘ > (x,6) , 

( 7 ) 
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and 



dw° 

dx 

dw° 

Rd& 


f 


(9) 

which describe the displacement field in general according to the Kirchhoff assumption, 
and eqs. (173a), (173b), and (166) which are the solutions for u°(x), v°(x), and w°(x) 
for 0 * N> N * , eqs. (174a), (174b), and (167) which are the solutions for u °(x ) , v°(x) , 
and w‘(x) for N = N * , and eqs. (175a), (175b), and (168) which are solutions for u °(xr) , 
v°(x) , and w°(x) for N<N* . Additionally, the assumption of Donnell was incorporated 
in the analyses of the previous chapters, namely, that for thin shells, the variable r can 
be replaced by the mean radius R. 

Under the conditions imposed by the Kirchhoff and Donnell assumptions, the partial 
differentials which make up the axisymmetric displacement gradients and rotations listed 


in eq. (211) simplify to 


du° d\>° 
t._ = — — z- 


dx 

2e = — 
' e dx ’ 


dx 1 ’ 


2r ze = = ~2r a , 




dx 


-2^1 = -2r 
dx 


2e xr = 0 , 






( 212 ) 


w 

e '* = T’ 


2 ‘e, - ~zr . 


e rr =0. 
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Note that under these assumptions, 


r xB ~ ~ e x6 » ^ r 9r ~ e 9r * 


(213) 


E. Specialization of the Three-Dimensional Equilibrium Equations and Boundary 
Conditions Under the Assumption of Axisymmetry 

The analyses previous to this chapter have been made under the assumption of 
axisymmetry. Since the solutions obtained in the previous chapters are to be 
incorporated into the three-dimensional equations derived previously in this chapter, the 
last form of the equilibrium equations and boundary conditions, i.e., eq. (208) and (209), 
will be specialized under the assumption of axisymmetry. 

Enforcing the conditions for axisymmetry, eq. (210), into the equilibrium equations of 

the previous section, i.e., eqs. (200), results in 


r|[(l * V"] 

- |/[(i * r ,y] - ° • 


(214a) 


r | lv’” + ( 1 ♦**)**] 


- r xr X 


x« 


+ 


ee 


+ t #r = 0 , 


(214b) 
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and 


* + ( 2 ^)' 8 ' * °1 
- 2 - (1 *e^c" - 2' t y = 0 . 


(214c) 


Enforcing eq. (210) on the natural boundary condition equations, eqs. (201), results 

[(1 + e„)o** + r xr t xr ]n x + [(1 + r ir t e ']n 6 

+ K 1 + «„) T ' r + r ,r a "]"r = i* S 2 (215a) 

OR u is specifiedon S 2 , 


and 


[( 2 **K + (1 + *ee) T * 8 ] n , 

+ [( 2 «xe)^ e + (* + «ee)o ee ]»e 
+ [( 2 «*e) TXr + ( ! + «ee) Ter ]»r = * # <® S 2 

OR v is specified on S 2 , 


(215b) 


and 


* (2‘e^ 6 * *”]«, 

+ [(- r „) T ” * (2«, r )T ,r * ®"]», - <' on S 2 

OR w is specified on S 2 . 


(215c) 


F. Simplification of the Equilibrium Equations and Boundary Conditions Through 
Elimination of Terms of Relatively Small Magnitude 

At this point, the three-dimensional equilibrium equations and boundary conditions 
have been simplified so that the solutions of the previous chapters can be used to evaluate 
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the displacement gradients, rotations, and intralaminar stresses at any point within the 
cylinder wall. It should be noted that these quantities are obtained from CLT. The 
solutions of the previous chapter are being used in order to later derive recursive closed- 
form relations for the interlaminar stresses at a given point along the length of the 
cylinder. 

In an attempt to further simplify the equilibrium equations, eqs. (214), and boundary 
condition relations, eqs. (215), a study will be conducted of the relative magnitudes of 
the terms which are being differentiated in the equilibrium equations. This is similar to 
the method implemented by Novozhilov (Ref. 6) to simplify the three-dimensional 
equilibrium equations based on the relative magnitude of the displacement gradients and 
rotations appearing in these equations in order to rationally obtain equations relevant to 
small deformations and small rotations. To be conservative, no assumptions will be 
made as to whether a term can be eliminated under a general case of cylinder stacking 
sequence, constitutive properties, geometry, boundary condition, or load level. 
Otherwise, the study will be conducted based on these parameters as they will occur in 
the subsequent calculation of the interlaminar stresses. This precaution is being taken 
due to the nonlinear nature of the CLT solutions presented in the previous chapters. 

The cylinders to be analyzed in the rest of this work are the same cylinders which 
were analyzed in the previous chapter on intralaminar stresses. Likewise, only clamped 
boundary conditions will be investigated. Thermally-induced preloading effects with 
AT=-280°F will be used in the CLT solutions. 

The maximum absolute values of the displacement gradients and rotations for the 
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cylinders subjected to compressive axial loads of N = 10% V*, N = 90% N * , and 
N - 99% N* are presented in Table III. As can be seen, the displacement gradients and 
rotations are considerably smaller than unity. 

Table HI. Maximum Values of Displacement Gradients and Rotations 


Cylinder 
& Load 

1 1 mix 

r=r mner 

1 ^xx 1 max 
^ ^ cuter 


■a 

IS 

1 1 max 

1 1 max 

[±45/0J2s 

N=10%N* 

.20*10- 3 

.29*10- 3 

.12*10- 2 

.11*10- 2 

0 

0 

[±45/0J 4T 

N=10%N* 

.20* 10- 3 

.31*10' 3 

.11*10- 2 

.n*io - 2 

.40*10" 5 

.10*10-* 

[(y T 45] 4T 
N=10%N* 

.27*10" 3 

.31 *10* 3 

.12*10" 2 

oo 

* 

.28*10" 3 

.17*10* 

[145/OJ* 

N=90%N’ 

.29* 10‘ 2 

.32*10" 2 

.12*10" 2 

.15*10-' 

0 

0 

[±45/0J 4T 

N=90%N* 

.29*10- 2 


.11*10- 2 

.14*10-' 

.28*10* 

.14*10" 5 

[0 2 /T45] 4t 

N=90%N* 

.35*10- 2 

.38*10" 3 

.13*10" 2 

.18*10-' 

.35*10^ 

.18*10" 5 

j| [±45/02l2S 
N=99%N* 

.37*10" 2 

.35*10- 2 

.15*10- 2 

.21*10-' 

0 

0 

[±45/0J 4T 

N=99%N’ 

.36*10- 2 

.32*10" 2 

.15*10- 2 

.20*10-' 

.30*10* 

. 18*10" 5 

[(y T 45 ] 4t 
N=99%N* 

.44*10- 2 

.42*10- 2 

.21*10- 2 

.25*10-' 

.39*10* 

.24* 10 s 


L Simplification of the First Equilibrium Equation 

The derivative with respect to x in the first equilibrium equation, eq. (214a), is 
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(216) 


|-[(1 ♦ <„)o“ * r„t»] . 

For the range of load levels investigated, the maximum absolute value of the 
displacement gradient e„ is less than Id 2 . Therefore, is small enough relative to 
1 such that it can be neglected, and the coefficient of o“ becomes 1. Since the 
magnitude of x xr is expected to be small relative to the magnitude of a xx , and since the 
rotation r XT is small relative to 1 for these load cases, it can be reasonably assumed that 
a xx is sufficiently large compared to r xr x xr to warrant the elimination of the latter term 

from the derivative with respect to x. 

The derivative with respect to r in the first equilibrium equation, eq. (214a), is 


+ r » a 1 • 


017) 


As with the derivative with respect to x, e xx can be neglected such that the coefficient 
of x XT is 1. Also, the interlaminar stress component a rr is expected to be small relative 
to x xr . Therefore, it can reasonably be assumed that x XT is sufficiently large compared 
to r xr a TT to warrant the elimination of the latter from the derivative with respect to r. 
In light of these results, the first equilibrium equation becomes 

eg + afr*") . o , < 218 > 

Bx dr 

which is a linear partial differential equation. 

2 Simplification of the Sec ond Equilibrium Equation 

The derivative with respect to x in the second equation, eq. (214b), is 
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(219) 


M 1+ • 

Referring to Table III, the maximum absolute value of e w for the load levels investigated 
is on the order of lCf 3 . Therefore, it can be neglected in the coefficient of x* 6 in the 
derivative with respect to x and the coefficient of x® r in the derivative with respect to 
r. The maximum absolute value of e x9 is less than 10 4 . Therefore, in the derivative 
with respect to x, the maximum absolute value of the coefficient of a xx is less than 10 4 
and the coefficient of x* e is 1. The intralaminar stress component x* e is, in general, 
one or two orders of magnitude smaller than a xx . Therefore, it can reasonably be 
assumed that x* e is sufficiently large relative to e x9 a xx in the derivative with respect to 
x such that the latter term can be eliminated. 

The derivative with respect to r in the second equilibrium equation is 

+ (1 + . ( 220 ) 

Based on the observations for the derivatives with respect to x, the coefficient of x 9r in 
the derivative with respect to r is 1 and the coefficient of x sr has a maximum absolute 
value of 10 4 . The interlaminar stress components x XT and x 9r are expected to differ by 
several orders of magnitude. Therefore, it is reasonable to assume that x 6r is sufficiently 
large relative to e x0 x xr in the derivative with respect to r such that the latter term can 
be eliminated. 

The terms in the second equilibrium equation that are not differentiated are 
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(221) 


Since the magnitudes of the intralaminar stress components o ee and x x6 are about the 
same, and the maximum absolute value of r XT is approximately 0.025 for the load cases 
investigated while the maximum absolute value of e 0r is less than 10 5 , it can reasonably 
be assumed that the term r XT x x ® is sufficiently large compared to e Br a w such that the 
latter term can be eliminated. Since the interlaminar stress component x 6r is expected 
to be several orders of magnitude smaller than x xe , the magnitude of the term r„x xe 

may be comparable to the magnitude of x . 

Based on these assumptions, the second equilibrium equation becomes 




( 222 ) 


which is a nonlinear partial differential equation. 



The derivative with respect to x in the third equilibrium equation, eq. (214c), is 


4 4 *1 • <223) 

The maximum absolute value of the displacement gradient e 9r is less than 10" and the 
maximum absolute value of the rotation r xr is 0.025. Therefore, the magnitude of the 
^ — r x ^j <r* x term can be assumed to be sufficiently large compared to the (2e I0 )x term 
such that the latter term can be eliminated. However, since x xr is expected to be several 
orders of magnitude smaller than o~, the magnitude of (~r xr )o xx can be expected to be 
comparable to the magnitude of x xr . Therefore, based on these observations, the 
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derivative with respect to x becomes 


^[(-r„)o“ ♦ x"] . (224) 

The derivative with respect to r in the third equilibrium equation is 

I'K-VK' + P*.,)*" * «"] • 025) 

Since the interlaminar stress components x xr and x 0r are expected to differ by several 
orders of magnitude, and in light of the magnitudes of the coefficients of these 
components, it can be assumed that (-r Xf ) x xr is sufficiently large compared to (2e 6r )T® r 
such that the latter term can be eliminated. The interlaminar stress component o rr is 
expected to be small relative to x xr . Therefore, the magnitude of the ( - r xr } x xr term can 
be assumed to be comparable to the magnitude of the a" term. Based on these 
observations, the derivative with respect to r becomes 

The terms in the third equilibrium equation that are not differentiated are 

• < 227 ) 

Due to the magnitudes of the absolute values of the displacement gradients e x6 and e 6r , 
the terms -2e xe x x6 and -2e 0r x Br can be neglected with respect to the term 
(1 + e ee j o®® . Also, due to the magnitude of e M relative to 1, it can be neglected in the 
latter term. 

Based on the observations of the relative magnitudes of the terms in the third 
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equilibrium equation, it can be simplified to 


r i * *1 + i r iw * °1 - = 0 • (228) 

d Simplification of the Natural Bou ndary Conditions 

Since the terms which appear as coefficients of the surface normals n x and n r in the 
natural boundary condition relations, eqs. (215), are identical to the terms being 
differentiated with respect to x and r, respectively, in the equilibrium equations, eqs. 
(214), the results of the previous sections where the equilibrium equations were 
simplified can be used directly to simplify the natural boundary condition relations. 
Therefore, based on the magnitude study for the equilibrium equations, the natural 
boundary condition relations become 

M*. * [(» * ««K* * v*1-. * M 5 ' - *' “ s > (229a) 

OR u is specified on S 2 , 

^ + [( 2 «xe) T * 9 + f 1 + *ee)° e6 ]"e + [* ir ]»r = ** 011 S 2 (229b) 

OR v is specified on S 2 , 


and 


R-U*- 


* *"]», 

* (2 «„)«•• ♦ t ,r ]5, 

+ o rr j n r = i r on s 2 


OR w is specified on S 2 . 


(229c) 


However, the coefficients of the surface normal n e must be analyzed as to the relative 
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magnitudes of the terms in these coefficients. In the coefficient of n e in the first natural 
boundary condition relation, eq. (229a), e xx can be neglected relative to 1 in the 
coefficient of r* e . Since the interlaminar stress component t 0 ' is expected to be small 
relative to the intralaminar stress component r* 0 , the term t* 6 can be assumed 
sufficiently large compared to r xr x 9r such that the later term can be eliminated. Based 
on these observations, the first natural boundary condition relation becomes 

[°“E + b x *K + b*'E = * x s 2 

(230a) 

OR u is specified on S 2 . 

In the coefficient of n 6 in the second natural boundary condition relation, is small 
relative to 1 such that it can be ignored in the coefficient of o ee . Also, since the 
maximum absolute value of e x0 is less than 10" 4 , it can be assumed that the term o 00 is 
sufficiently large compared to the term (2e ie )T* e such that the latter term can be 
eliminated. Therefore, the second natural boundary condition becomes 

+ [o®°] w e + = * 8 on S 2 

(230b) 

OR v is specified on S 2 . 

In the coefficient of n 6 in the third boundary condition relation, the term 
can assumed to be small relative to the other two terms, since e 6r is small. However, 
although t 0r is expected to be small relative to the intralaminar stress component t* e , 
and since the maximum absolute value of r xr is approximately 0.025, the magnitude of 
the term (-r^r* 0 can be assumed comparable to the magnitude of the term x 0, ‘. 
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Therefore, the third natural boundary condition relation becomes 


K,) 4 *1". 

* 4 

+ [(' r „) T " + «"]», * >' 0,1 S 2 

OR w is specified on S 2 . 


(230c) 


In summary, after simplification by the consideration of the relative magnitudes of 
the terms appearing in the equilibrium equations under the assumptions of axisymmetric 
response and loading, and the assumptions of Kirchhoff and Donnell, these equations are 

r da xx + d(r* xr ) = 0 > (231a) 

ctx dr 


+ . r ^> * . o , 

cbt dr 


(231b) 


and 


4 1 -] * 4 - ° 


(231c) 


By the same method, the natural boundary conditions become 


[<r**]n x + [t x8 ]n e + [T xr ]», = i x on S 2 

OR u is specified on S 2 , 


(232a) 


and 


ff]*. 4 l""]*. 4 = *• “ S 2 

OR v is specified on S 2 , 


(232b) 


143 


and 


+ [(- r *rK + ^]»e 

+ [(~ r xr)^ Xr + »"]», " *' ° nS 2 

OR w is specified cm S 2 


(232c) 


In the following chapter, these equations will be solved and the interlaminar stresses 
will be calculated for the three cylinders of interest. 
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VL SOLUTION OF THE THREE-DIMENSIONAL EQUILIBRIUM 
EQUATIONS FOR THE INTERLAMINAR STRESSES 


In chapter five, the three-dimensional equilibrium equations and boundary conditions 
were derived for a linear elastic body in cylindrical coordinates for finite strains, under 
the assumption of axisymmetric loading and response, and under Kirchhoff s and 
Donnell’s assumptions for the displacement variables. These steps were performed in 
order to obtain three-dimensional equilibrium equations and boundary conditions 
compatible with the assumptions and solutions of the previous chapters, which concluded 
with the derivation of the intralaminar stress relations. Under these conditions, the 
closed-form intralaminar stress relations can be used in the solution of the interlaminar 
stress components, analogous to the method implemented by Pagano (Ref. 3). 

In this chapter, the first of the equilibrium equations, i.e. 

do" + d(rx xr ) = Q # (233) 

dx dr 

will be solved for the interlaminar stress component z xr . Since o xx {x,f) is known at any 
point along the length of the cylinder and at any point through its thickness from closed- 
form relations obtained previously, the partial derivative with respect to x of a xx is also 
known through analytical differentiation of the closed-form solution to o"(x,r) . 

The interlaminar stress components x 6r and o" are assumed to be small relative to 
the interlaminar stress component t xr . For an axisymmetric response, t would be 
expected to be small. Other researchers (ref. 7) have found o" to be small for 
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the case of cylinder bending, and that is assumed to be the case for the cylinder 
compression problem studied here. Therefore, t 8r and a rr will not be solved for or 
calculated in this work. However, the solution method to be presented for the 
determination of the interlaminar stress component t xr from the first equilibrium 
equation is directly applicable to the determination of the other two interlaminar stress 
components. 


A. Solution of the First Equilibrium Equation 

The differentiation implicit in eq. (233) can be distributed, and the result simplified 
to give 


dx*' + 
dr r 


da xx 

dx 


( 234 ) 


As mentioned previously, the partial derivative appearing on the right-hand side of this 
equation can be calculated from the closed-form relation for a xx (x,r ) . For a given layer, 
the intralaminar stress component o xx (jc,r) is given by eqs. (19) and (20), which are the 
CLT stress-strain relations, and eqs. (12) and (128), which are the axisymmetric 
kinematic relations. They are repeated here for convenience: 

CLT Stress-Strain Relations: 


°x = Qn e * + Qi2 c e + “ °i 

°a = Ql2 C x + Q22 C e + ^26**6 “ °$ 

X *e = Ql6 e jr + + QmYxO T j« » 


where 


( 19 ) 
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( 20 ) 


o T x = (Qn«, + Qi2 a e + Qi6 a ^) Ar 
0g = (Qn«, + Q22 a e + ^26*^)^^* 

tJe = (Qi6«x + Q26«e + Q<* a *)* T - 


Axisymmetric Kine matic Relations; 

0 0 

E I = t, + Z>=i 

0 0 

Cg = Cg + Z K® 

Y* = Y* + Z <e , 

(12) 


where 


PS=« 


*: + ip 

dx 2 




. *p: 


1 o W° o dv° 
kJ»0;<,=0. 


(128) 


In order to proceed with the solution process, the first of eqs. (19) must be 
differentiated with respect to x. Since the transformed reduced stiffnesses C? i; are 
assumed constant along the length of the cyUnder, the partial derivative of a xx with 


respect to x, in terms of the strain components is 


do” _ r, 8e ' + O ♦ O ^ (235) 

1 iT ' Q " aT °»1T at 


Therefore, the partial derivatives of the inplane strain components given in eqs. (12) and 
(128) are required. Performing the differentiation results in 
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dw d 2 ^ 
dx dx 2 


d 3 w° 



d 2 u° 

dx 

dx 2 

de 0 

1 dw 

dx 

R dx 

a Y, e 

_ d 2 v° 

dx 

dx 2 



(236) 


These partial derivatives involve the derivatives of the solutions for the reference surface 
displacements u°(x), v°(x), and w°(x). For instance, the first of eqs. (236) involves 
the second derivative of u °(x ) . This term can most easily be derived from the definition 
of N x for the axisymmetric problem, given in eq. (134) as 

", - A u< * - N, T . (237) 

where t° x is given by eq. (13) as 


du° , 1 a o 2 du° . 1 dw° 
e, = + — p r = + — 

dx 2 dx 2 1 dx 


\2 


(238) 


du° 

Substituting eq. (238) into eq. (237), and solving for results in 

dx 

= _L(w . if **:] 2 

dr A ' x * I A. n R A -»-* “)l I 


A n dx 2 2( dx ) 


(239) 


Differentiating the above expression once results in 


_ -^12 dw° + ^ii d*w^_ _ dw° dhv^ (240) 

dx 1 A n R dx A n dx 3 dx dx 2 
Since the solution for w°(x) is known and it is continuous, its derivatives are obtainable 
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and the expression in eq. (240) is known for any point along the half-length of the 
cylinder. 

dt x . 

Substituting eq. (240) into eq. (236) results in an expression for — in terms of 
derivatives of w°(x), i.e., 


H 

U 

^12 dw° 

, *n d*w° 

dw° d 2 w° 

dw° 

+ 

d2 ”\ z 


dx 

i4 u /? dx 

A n dx * 

dx dx 2 

dx 

dx 2 

, dx* j 

dt x _ 

A l2 dw° 

, B n d*w° 





dx 

i4 u /J dx 

A n dx * 

1 dx* J 





(241) 

Since eq. (234) involves differentiation with respect to r, it will be necessary to express 
the above equation in terms of r instead of z. This is easily accomplished by substituting 
the definition of the local z coordinate, i.e., 

«-r-*, <242) 

into the last expression in eq. (241), which results in 

= ±2 dw° Jjn +R )d^ +r (^l\ (243) 

dx A n R dx [i4 n ) dx * { dx * ) 

The partial derivative of in eqs. (236) involves the constant R , the mean radius 
of the cylinder, and the second derivative of w°(x), which can be obtained through 
differentiation of the solution for w °(x) , derived in chapter 3. 

The partial derivative of Y x e in eq* (236), i.e., , can be obtained by 

differentiating the equation for y x9 in terms of w®(x) given by eq. (135b), namely, 
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(135b) 


y° - ( d\> 0 J_) 

* { A <* dx 2 A" 

Differentiating this equation once with respect to x results in 

dv<e _ ^i6 d 3 w* 

& x A 66 dx 3 

Substituting the expressions for the partial derivatives of the inplane strain 
components with respect to x , eqs. (243), (244), and (236), into the expression 
describing the partial derivative of a xx , eq. (235), results in 


do" _ 

r 

O 

A 12 

dw° 

j . 


d 3 w° 

dx 

^11 

A U R 

dx 

Un J 

dx 3 


1 dw° 

+ o 

^16 d 3 w° 


R dx 


A 66 dx 3 




' * 
"0 


“<?u 

1 



This expression can be rewritten in the form 


do xx _ 

[<*.“*] 

► 4. » 4 


dx 

| dx 

T r 1 

[ dx 


where 


do xx “ 

dx 



dw° 

dx 





1 dw° 

+ n 

^16 d 3 w° 


R dx 

+ V W 

* 



(245) 


(246) 


(247a) 
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and 


(247b) 


do 3 


dx 


•\-Qii 


d 3 *' 


dx r 3 


Therefore, at any axial location along the length of the cylinder, the derivatives defined 
in eqs. (247) can be calculated from the derivatives of the solution for w°(x) and the 
transformed reduced stiffnesses Q tj for the layer which corresponds to the coordinate r . 
It should be noted that the derivatives in eqs. (247) are known functions of the x 
coordinate. In order to make this point clear and to simplify the notation, the definitions 


m(x) = 


do « M 1 

n 

A 12 

dw° 

(B i, ) 

— +R 

d 3 w 0 ' 

dx 

Vll 

1 

»— » 

X 

dx 

^ii j 

dx 3 


1 dw° 

+ n 


d 3 w°' 


R dx 

Vie 

i — 

dx 3 



(248a) 


and 


m 


do*^_ 

dx 



(248b) 


are introduced. Substituting eqs. (246) and (247) into eq. (234) using the definitions in 
eqs. (248) results in the partial differential equation 


dt xr (xs) + **'(* £) = _ (m(x) + r[fc(x)]} . @49) 

dr r 

For a specified temperature change, boundary condition, axial load level N x , and axial 
position x=x, the terms m(x) and b(x) are known quantities. Therefore, at a specified 
axial location x =x, eq. (249) becomes a nonhomogeneous ordinary differential equation. 

That is, 
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(250) 


</T* r (s,r) + x >r (jc > r) 
</r r 


= -{ro(x) + r[&(*)]} 


Equation (250) is of the form 


dF F . 

— + - = 5 + hr , 
dr r 


(251) 


where the substitutions 


F = Fir) = T* r (*,r) ; $ = -mix ) ; h = -i(J) 


(252) 


have been made. By introducing a change of variable according to 


r * 


(253) 


eq. (251) becomes 


4 ^ + F = ge z + he 2t . 


(254) 


This differential equation has a homogeneous solution of the form 


or F_ = y . (255) 

where A is an unknown constant to be determined, and the particular solution 

F part = set + 2tg2Z * (256) 

Substituting eq. (256) into eq. (254) results in 

(25)e l + (3f)tf 2l =^e l + Ac 2t - 2s = g ; 3f = h . < 257 ) 

Therefore, the particular solution is 
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or 


(258) 


F = $-e z + — e 2 * 
^ 2 3 


F w = — r + — r* 

port. 2 3 


Hence, the solution to the nonhomogeneous ordinary differential equation (251) is 

F(r) = i + ir + V . (259) 

w r 2 3 

Since the step relating the partial differential equation (249) to the ordinary differential 
equations (250) and (251) was taken by restricting eqs. (250) and (251) to a particular 
x location, it must be realized that the constants A, g, and h appearing in eq. (259) are 
unique for each x location. Therefore, g and h in eq. (259) vary along the length of the 
cylinder and are known quantities obtained by the relations in eqs. (248), and A -A(x) 
is a unique constant for each x location, which will be determined later. Thus, by using 
the definitions from eq. (252), and noting that A =A(x) , for a particular axial location 
x = x, the solution presented above in eq. (259) becomes 

t xr ( - r) = Mil _ ”»(*) r _ kill r 2 . (260) 

r 2 3 

It should be noted that since the stress-strain relations of eq. (19) apply only within a 
given layer, the solution given in eq. (260) also applies only within a given layer. To 
make this distinction clear, a superscript (k), which denotes the layer number from 
k = l at the inner layer to k = K at the outer layer, will be used from this point onward. 
For example, eq. (260) will now be written 
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9 


(261) 


T *r»(- r ) = A W (x) _ IW < * ) (J) r _ ^ W (x) r 2 
r 2 3 

which is valid in the range 

r (*-i) s r $ r (*) 


(262) 


Therefore, for a cylinder with K layers, there are K constants A (k) (x) to be determined 
in as many equations (261). In eq. (262), the superscript on r can take on the range 
from 0 to K, with r m = - yj being the radius of the inner surface of the cylinder, 

and r w = |/? + j being the radius of the outer surface of the cylinder. Therefore, r (k) 
corresponds to the k A layer interface. 

Thus far, the pertinent kinematic relations, constitutive relations, and the first 
equilibrium equation have been used to derive the solution form of t xr . Next, the 
boundary condition relations of eqs. (232) must be satisfied for each layer of the 
cylinder. However, this is a trivial matter since the displacements are specified 
everywhere within the cylinder and on its boundaries due to the assumptions of the CLT 
analysis. 

In this analysis, the adjacent layers are assumed to have a perfect bond at the layer 
interface. This implies that all of the displacement components and interlaminar stress 
components are continuous across the interfaces of adjacent layers. As noted in the 
previous paragraph, the displacements are continuous across the interfaces since they 
have been prescribed to vary linearly through the thickness of the cylinder wall in the 
CLT analysis of chapters 2, 3, and 4. However, the interface condition for the 
continuity of the interlaminar stress components has not yet been addressed. In order to 


154 



solve for the unknown constants A (t) (J), the interface continuity of the interlaminar 
stress component t ir will be used. For the *=1 through K layers, there are (tf-1) 
interface continuity conditions and two surfaces on which x Jtr can be specified, namely 
the surfaces at the inner and outer radii of the cylinder. This results in (K + 1) conditions 
from which the K unknown constants A ik \x) can be determined. Since the interface 
conditions for x xr must be satisfied in order to comply with the assumption of perfectly 
bonded layers mentioned previously, the boundary condition for x xr at either the inner 
or outer radius will have to be ignored. Since the axial loading investigated thus far in 
this work does not consist of an applied traction on either the inner or outer surfaces of 
the cylinder, the condition 

- 0 , (263) 

will be imposed. Therefore, for a particular axial location, x=x, the equation 
describing x ZT for the first layer at the inner radius can be written as 

T ,r°>/r ,w\ . ^ (1) (*> _ £^ m O) ( j) _ &-b m (x) = 0 . (264) 

v * ' r (0) 2 3 

This equation can be solved for j4 (1) (x ) to give 

. < 2<S > 

2 3 

where the constants m (,> (3c) and b (t) (x ) are calculated from eqs. (248) using values for 
the transformed reduced stiffnesses Q i} for the first (inner) layer of the cylinder. 

The interface continuity condition for x XT is 
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(266) 




(x) = T"*(J) 


over the range of k = 1 to K-l . Substituting eq. (261) into eq. (266) results in 


A^\x) 

r W 


r<*) 


\x) 


■m 


A^(x) 

r (t) 


r (t) 

~2 


m (k \x) - kf£b<*)(Z) 


9 


(267) 


which can be rearranged to provide a recursive relation for A (k * l) (x), i.e., 


A (t+1) (x) = 

A (k) (x) + k^lm^ix) - m<*>(x)] + - b (k \x )] 

2 3 


(268) 


Hence, the K constants A (k) (x) in the equations describing T* r (x J) are known quantities 
and the interlaminar stress component x xr can be calculated at any radial location r at 
a specified axial location x=x where the CLT solutions have been previously calculated. 


B. Calculation of the Interlaminar Stress Component r" 

In the following, numerical results will be presented which illustrate the relations 
between the interlaminar stress component T Jtr and the radial coordinate r for the three 
cylinders analyzed in this investigation. Clamped boundary conditions are enforced, 
compressive axial loads of N = 10%//* and N = 90% N* are applied, and thermally- 
induced preloading effects corresponding to a temperature change of -280°F are included 
in the CLT solutions used for the calculation of the terms of the solution of x xr involving 
the intralaminar stress component a xx . As noted in the previous section, the solution for 
t (t) (jc,r) is calculable only after an axial position x=x has been selected and the 
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derivatives appearing in the definitions of the "constants" A ®(x ) , wi ( *(■* ) » snd b ( x ) 
have been calculated for each layer. The axial positions investigated will be: (a) The end 
of the cylinder at x=+L/2; (b) The axial position x where the rotation r XT has a 
maximum absolute value, and; (c) The axial position x where the interlaminar shear 
stress resultant, Q x , has a maximum absolute value. An expression relating Q x to the 
CLT solution is contained in the axisymmetric version of eq. (126b), namely, 


dM, 


Q = — i + ™x~T~ 

Vjt dx x ^ 


+ N. 


dw 


(269) 


dx 


The term ^ can be calculated from the expression for M x in eq. (134), i.e., 
dx 


u, - 8„«; * ♦ B„rU - D u ^ - Ml , 


2 ,,,o 


dw 


(270) 


dx 2 


which, by substituting the expressions relating t° x and Y,e to solution for w°(x) and 
its derivatives, i.e., eqs. (135), and differentiating the result once with respect to x 

becomes 


dM x 

B n *16 _ D 

d 3 w° A 

f B n B n A l2 ^ 

dw° (271) 

dx 

— + — - i/ii 

A A 11 

[ A l\ ^66 J 

dx 3 

R A„R j 

dx 


The linear form of Q x is obtained by eliminating the second term in eq. (269), i.e., 



(272) 


as given by equation (271). Since the solution for x xr was obtained from the linearized 
form of the first equilibrium equation, then it is expected that the formal definition of the 


resultant Q x , i.e, 
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should agree with the expression for Q x given in eqs. (271) and (272). 
Therefore, the axial position x corresponding to the largest absolute value of 



_ii + -D 

d 3 w° 

-L 

' B\2 ^h-^12 

dw ° (274) 

A A 11 

"66 ) 

dx 3 

R A n R J 

dx ’ 


will be one of the axial locations used to calculate the relation between x xr and r. 
Hence, as a check of the interlaminar stress calculation for t xr , the values of x xr will 
be numerically integrated along the radial direction as shown in eq. (273), and the result 
will be compared to that calculated through eq. (274). 


C. Numerical Results for the Interlaminar Stress Component T xr ; Case of 
Thermally-Induced Preloading Effects and a Compressive Axial Load with 
Clamped Boundary Conditions 

Fig. 41 through Fig. 45 illustrate the relationship between the normalized interlaminar 
shear stress and the radial coordinate p for the compressive axial load levels 
N = 10% N* and N - 90% N* at various locations along the length of the three cylinders. 
The shear stress x xr has been normalized by the quantity (NfH ) , as the intralaminar 
stress components were in chapter 4. The radial coordinate has been redefined as p , 


where 



( 275 ) 


r - R _ Z_ 
H H 


and 


-0.5 s p s +0.5 . 


(276) 


Each figure represents the interlaminar shear stress response of the three cylinders for 
the two load cases and a particular axial location. In each of the figures, horizontal grid 
lines and symbols are used to designate the 15 layer interfaces. Fig. 41 illustrates the 
shear stress response at the end of the cylinders, x=+L/2, associated with the low axial 
load level. While distribution of the shear stress in the symmetric laminate is shown to 
be symmetric about the mean radius, p =0, the distributions of the shear stress in the 
two unsymmetric cylinders are skewed, with the shear stress having larger magnitudes 
to either side of the mean radius. For instance, the [(y-45/ +45] 4T laminate has a peak 
shear stress magnitude at a radius two layer thicknesses inside of the mean radius. This 
behavior leads to larger slope discontinuities in the interlaminar shear stress at the 
interface of layers of differing orientations on the side of the mean radius to which the 
response is skewed. Conversely, the slope discontinuities are smaller to the other side 
of the maximum response location. Due to the opposite signs of the B matrix terms, the 
[+45/-45/0J4T laminate has a peak interlaminar shear stress at a radial location outside 
of the mean radius. These trends can be observed for the higher load level as well, i.e., 
the case represented in Fig. 42, which also represents the interlaminar stress component 
at the end of the cylinder. However, partly due to the scale used to illustrate the higher 
load level response plots, the effects appear to be smaller. It is also seen that the shear 
stress at the end of the cylinder for the [(V-45/+4514T laminate is larger for both the low 
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and high load levels then for the other two laminates. 

For the low load level, the axial location of the maximum shear stress resultant Q x 
is at the ends of all three cylinders. However, for the high load level the axial location 
of the largest Q x occurs at x/L=0.480 for the [+45/-45/0J2S and [O 2 / -45/ -F 45] 4T 
cylinders and at x/L=0.478 for the [+45/-45/02 Ji T cylinder. The interlaminar shear 
stress response for this case is presented in Fig. 43. Again, the [0 2 /-45/+45] 4T cylinder 
has the largest shear stress magnitude. Shown in this plot is the fact that the shear stress 
of all three cylinders at this axial location is also skewed relative to the mean radius of 
the cylinder. The reason for this behavior is not immediately apparent, but it could be 
a function of the magnitudes of the first and third derivatives of w °(x) at axial locations 
away from the end of the cylinders, as they appear in the relations for the solution for 
T* r , relative to the magnitudes of these terms at the ends of the cylinders. 

Fig. 44 and Fig. 45 represent the interlaminar shear stress response at axial locations 
for which the cylinders possess peak magnitudes of the rotation, r xr , for the low and 
high load levels, respectively. These axial locations are in the vicinity of x/L=.470 for 
all three laminates at the high load level, and at x/L= 0.480 for the [+45/-45/0J2S and 
[0 j/- 45/+45] 4T cylinders and at x/L=0.5 for the [+45/-45/()J 4T cylinder at the low load 
level. Again, it is noted that the response of all three cylinders is skewed relative to the 
mean radius. For the shear stress results previously discussed, i.e., at the low and high 
load level and where the axial locations were that of maximum Q x and the cylinders’ 
end, the [OjMS/ +45]^ cylinder has the largest overall shear stress magnitude of the two 
three cylinders. This is also true for the high load level where the axial location is that 
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of the maximum rotation r xr , as seen in Fig. 45. However, for the low load level and 
at the axial location of maximum rotation r xr , the overall magnitude of the shear stress 
for the [+45/-45/0j4 T cylinder is larger than the overall magnitude for the [0 2 /-45/+45] 4T 
cylinder. This reversal in trend indicates that the thermally-induced preloading effects 
dominate the response at the low load level. This characteristic was observed with the 

interlaminar stresses discussed in chapter 4. 

The responses shown in Fig. 41 through Fig. 45 where numerically integrated using 
Simpson’s 1/3 quadrature rule. Since the relation describing r xr as a function of r is 
parabolic in r, and Simpson’s quadrature involves approximation of the data by a series 
of parabolic segments, the results of the numerical integration by Simpson s rule are 
independent of the number of data points used. The results presented in the figures were 
numerically integrated on a layer-by-layer basis and the integration results were summed, 
as indicated in eq. (273). Also, the linearized relation for the shear resultant Q x , given 
by eq. (274) was calculated at the axial location corresponding to that at which the x 
data was obtained for each relation in each figure. These steps were conducted as a 
check of the accuracy of the solution for x xr . The results of these calculations are 
presented in Table IV. The right-most column of Table IV is the CLT relation for the 
shear resultant Q x as given in eq. (274), and the second column from the right is the 
result of Simpson’s quadrature on the t* r data. It is seen that there is excellent 
agreement between the two relations for all load cases and axial locations investigated. 
This result suggests that the solution for t" derived from the linearized version of the 
first equilibrium equation, eq. (233), is accurate for these cylinders under the load levels 
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and boundary conditions studied. 

This chapter brings to a close the investigation of the displacement and stress 
response of unsymmetrically laminated cylinders. Though specific cylinders were 
considered, general conclusions can be drawn from the results presented. A discussion 
of these conclusions is the subject of the final chapter. 



T* r / (N/H) 

Fig. 41. versus p at x/L=0.5 for [+45/-45/0J2S, [+45/-45/0J 4T , and 

[(V-45/+45] 4T Cylinders, N=10%N\ AT=-280°F. 
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r* r / (N/H) 

Fig. 44. versus p at x/L locations of maximum r„ for [ +45/-45/0 ,L <! , 

[+45/-45/0J4T, and [O 2 /- 45 / -t-45] 4T Cylinders, N=10%N\ 
AT=-280°F. 
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T* r / (N/H) 

Fig. 45. t° versus p at x/L locations of maximum r„ for [+45/-45/0 2 ] 2s , 
[+45/-45/02]4 T , and [02/-45/+45]4 T Cylinders, N=90%N*, 
AT=-280°F. 
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Table IV. Comparison of the Shear Resultant Q x as Calculated from the Interlaminar 
Shear Stress r" and the Derivative of the CLT Relation for M x . 


Cylinder 

Stacking 

Sequence 

Axial Load 
Level, 
N/N* 

Axial 

Location 

Criterion 

Axial 

Location, 

x/L 

! x * r do 

dMx 1 

dx 

J (N/H) dp 


10% 


0.5 

0.0157 

0.0158 


10% 


0.5 

0.0149 

0.0149 

[0 2 /-45/ +45] 4T 

10% 


0.5 

0.0182 

0.0183 

SBSBSM 

10% 


0.477 

0.00572 

0.00572 


10% 

m 

0.5 

0.0149 

0.0149 

DRii 

10% 


0.485 

0.00909 

0.00913 

mmm 

90% 


0.5 

0.00524 

0.00524 


90% 


0.5 

0.00454 

0.00456 

[(V-45/+45U 

90% 


0.5 

0.00650 

0.00654 


90% 


0.470 

0.00846 

0.00851 


90% 


0.468 

0.00798 

0.00803 

[0,/-45/+45] 4T 

90% 

msm 

0.470 

0.01010 

0.01014 


90% 

KgSRii 

0.480 

0.00954 

0.00954 

HBfiHSPI 

90% 


0.478 

0.00888 

0.00892 

1 [0 2 /-45/+45] 4T 

90% 

■m 

0.480 

0.01119 

0.00127 


f i xr 

1 normalized to correlate with I dp 

J ( N/H ) 


165 


























































vn. CONCLUSIONS 


A. Discussion of Displacement Results 

In chapter 3, results were presented for the mid-plane displacements v°(jc), 

and w°(x) as they vary along the length of one symmetrically laminated cylinder and two 
unsymmetrically laminated cylinders with equal and opposite B matrix terms. From the 
results for the preloading response due to thermal effects, it was observed that all of the 
cylinders expand axially due to the negative value of the effective axial thermal expansion 
coefficient of the cylinders and the negative temperature change, AT, from consolidation 
temperature to ambient temperature. In addition, all three cylinders respond to the 
temperature change by a reduction in the radii of the cylinders. For the symmetric 
cylinder, the tangential displacement is zero and radial displacement is constant along the 
cylinder’s length. For the unsymmetric cylinders, the tangential and radial displacements 
vary along the cylinders’ length, particularly near the ends. These two cylinders’ ends 
"curl" radially and twist, producing a boundary layer effect. The direction of the curl 
for the cylinders is opposite, depending on the sign of the thermally induced moment 
A#/, which is a consequence of the opposite signs of the B matrix terms for these two 
unsymmetric cylinders. Any unsymmetric laminate will exhibit this characteristic. This 
thermally induced curl and twist must be accounted for in the boundary conditions for 
subsequent axial loading and could have an effect on the onset of buckling of the 
cylinders. This is because the radius at the end of an unsymmetric cylinder through 
which the axial load is applied is different from the radius away from the ends. 
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Displacement results of the case of compressive axial load, including the thermally- 
induced preloading effects discussed above, indicate that the thermal effects have a 
measurable influence on the radial displacements under axial load. Neglecting the 
thermally-induced preloading effects results in smaller radial deformations. For a highly 
unsymmetric stacking sequence, the effects become larger. The load-induced axial and 
tangential displacements appear to be relatively unaffected by the thermally-induced 
preloading effects. 

The shape and magnitude of the tangential and radial displacement responses change 
as the level of the compressive axial load increases. In particular, the length of the 
boundary layer, where the variation of the tangential and radial displacements fluctuate, 
increases as the load level increases. The axial displacement remains virtually linear 
along the axial direction, even at the high load level. 

Comparison of the displacement results for simply supported and clamped boundary 
conditions reveal that the simply supported case yields a larger range of tangential and 
radial displacements for each cylinder although the axial displacement response is about 
the same. 

An important conclusion is that if unsymmetrically laminated cylinders are to be 
analyzed or manufactured, it is important to include the thermally-induced preloading 
effects not only in the prediction of the overall displacement behavior of the cylinder 
under compressive axial loading, but also in the specification of the end conditions for 
both the analysis and the fixture design. The extent to which these effects are important 
depends on the material properties, stacking sequence, and the axial load level. 
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B. Discussion of Intralaminar Stress Results 


An interesting observation was made of the intralaminar stress results in chapter 4. 
At the low load level investigated, N = 10% N* , the fiber-direction stresses, o u , in the 
±45° layers is larger than the fiber-direction stresses in the 0° degree layers of the 
[+45/-45/0J2S. [+45/-45/0J 4T , and [(y-45/+45] 4T cylinders. This is surprising since 
the 0° layers have fibers aligned with the axial load and are expected to bear the majority 
of the axial load. This result indicates that at low load levels the thermally-induced 
preloading effects dominate the fiber-direction intralaminar stress response. The 
intralaminar stress component perpendicular to the fibers, o n , and the intralaminar shear 
stress component, t 12 , are small relative to the fiber-direction stresses. 

When the axial load level is increased to N = 90% N * , the thermally-induced 
preloading effects are seen to be in the background, based on the observation that now 
the 0° layers have a larger fiber-direction magnitude than the ±45° layers, as would be 
expected for a cylinder subject to large axial loads. The intralaminar stress component 
perpendicular to the fibers, , is observed to be compressive. This result virtually 
eliminates the potential for matrix micro-cracking due to this stress component. The 
possibility of inplane shear failure also seems low, since the magnitude of is 
observed to be low, even for this high load level. It should be noted that the application 
of the relations for the principal material stress components presented in chapter 4 to 
conventional CLT plane-stress failure theories can easily be accomplished, although it 
was not done as part of this work. 
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C. Discussion of Interlaminar Stress Results 

In order to more fully develop the relations for the material failure modes of the 
composite cylinders, a solution for the interlaminar shear stress x xr was developed and 
used to calculate the interlaminar shear stress for the cylinders and load cases 
investigated for the intralaminar stresses. It was reasoned that the remaining interlaminar 
stresses x* r and a" are small relative to x xr in the context of axisymmetric loading. 
The solution for the interlaminar shear stress x xr was derived through a rational 
simplification of the three-dimensional equilibrium equation for the axial direction. 
Based on the excellent comparison between the results of numerically integrating the 
interlaminar shear stresses and the relation for the linearized shear stress resultant Q x 
from the CLT solution of the third chapter, it is recognized that the solution obtained for 
the intr alam inar shear stress was accurate. It was hoped that this would be the case since 
the method used to simplify the three-dimensional equilibrium equation, and the solution 
used to obtain the other stress component appearing in the equilibrium equation, o a , 
both u tilize d the assumptions of Kirchhoff and Donnell. 

It is observed that at the ends of the cylinders, the interlaminar shear stress response 
of the symmetric cylinder is symmetric with respect to the mean radius, R, or p =0. 
The two unsymmetric cylinders, with [02/-45/+45]4 T and [+45/-45/0J4T stacking 
sequences and equal but opposite B matrix terms, have peak interlaminar shear stresses 
at radial locations to the inside and outside of the mean radius. However, the shape of 
the in terlam inar shear stress response through the cylinder wall is seen to vary along the 
length of the cylinder. In particular, for both the low and high load levels, the location 
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of peak interlaminar stress for the symmetric cylinder is seen to move away from the 
mean radius as the axial location at which the response is calculated is varied from the 
end of the cylinder, x/L=0.5, to the axial location of peak rotation, r xr . This is also 
true as the axial location is varied from the end of the cylinder to the axial location of 
peak Q x , for the high load level. 

It is also observed that the axial location of peak Q x occurs at the end of the 
cylinder, x/L=0.5, for the low load level, while the axial location of peak r xr is at 
x/L — 0.48 for the high load level. For this load level the interlaminar shear stress 
response is skewed relative to the mean radius. 

Comparing the results for the interlaminar shear stress calculated at the axial location 
of peak rotation, r xr , for the low and high load levels, reveals that the [+45/-45/0J 4T 
cylinder has the largest shear stress for the low load level, while the other unsymmetric 
cylinder has the largest shear stress for the high load level. This difference is yet again 
an indication of the thermally-induced preloading effects dominating the cylinder response 
at low load levels. 
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